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This  technical  documentary  report  has  been  reviewed  and  is  approved. 


HOLLAS?  B.  ITVM2S,  JB. 

acting  CMcf ,  Structures  Llv'sicr. 

IF  *: ig-t  Dyris'.cj  Laboratory 


a&3TR*CT 


The  equations  of  eotloc  applicable  to  the  Sir-Pi £Tfct -of -Credos 
Structura!  r- Program  (CLP)  are  derived  in  this  report.  These  equa¬ 
tion  written  for  the  determination  of  vehicle  structural  .oods  and 
response  due  to  aerodynamic  loads,  *oada  due  -*  control  junace  deflec¬ 
tions,  *na  <nviron*eatal  disturbances.  Arbitrary  elastic  degrees  of  frj^doa 
(ving  terdiug,  wing  torsion,  body  bending,  etc.)  and  rue)  cKsh  eruatioac 
are  inc<  .7 orated  Into  tne  overall  analysis. 

Newtonian  flov  theory  it  used  for  obtaining  idealised  s-'xiyrumlc 
pressure  distributions  since  it  is  the  simplest  aerodynamic  theory  that 
offers  sufficient  generality.  Accelerations,  deflections,  sbo&r  forces 
and  bc&iiag  *x>oeats  at  arbitrary  station*  can  be  computed. 

This  technical  docuaeotary  report  has  been  reviewed  and  is  approved. 


7.0UJCO  3.  USSX S,  JR. 

Acting  CM*?,  Structure*  Siv'sicr. 
\F  z'  ight  Dynasties  Laboratory 


1  I5C5TT«xmOH 


TABU  OF  COSTZXT^ 


£.  K2ASUT1C3  OF  72E  VZSICLE  A3D  BASIC  ASSUCPTICSS .  k 

2.  7CGC/2S.  HCICSftS,  AW)  THOWdC  RALAMCIlC" .  10 

s*  •JkUATicFS  o?  Mane*  rc«  tzk  a^snc  detcwocicks . ic 

5.  vcaa c  ajd  aaajT  khjoiok . 2J* 

6.  nucncAL  extpxssig*  or  the  ikrtial  rcfwuus .  26 

7.  TACTICAL  EOTSSSlC*  07  W5  AJfcODnAKlC  FOfcWLAS .  kl 

0.  IfCLOBlO*  or  THFUST  i CROSS  . .  u 

9.  BBUCTICK  (XSIMC3  07  ICWaSLS  STRUCTURAL  SBCTIOSi .  k< 

10.  7CPKULAS  TOR  B!2  SBPJOTJRAL  IOADS  .  k7 

11.  XJCLU3X0I  07  7USL  SI06HMG  K77SCT3 .  J*9 

12.  POIM  07  ROTAXICS?  hS D  EfaWQC  3AL«CIS<5  07  CSKCAIK  DATA  TOR 

TUHJOJE  SKTTlCWiJ  .  5? 

13.  ROTunrcss .  59 

APPZXDXX  I  -  BASIC  KBSfJLASTOIS  TDK  ?JSL  SUiSHIJC .  60 

AFFSXDU  11  -  STWOLS,  I^TA  TO  &S  SUBKTVi3K),  COOWfiXOaS  AJTO 

SVimOCa  OSKD  IK  IKS  STSUCTORAL  LOAUJ  ITOGRaM  .  .  &) 


W 


ILLUSTRATIONS 


*.  vfordliwkte  3y*tcaa . ^ 

2.  * cLtcie  tad  JWtioe  Coordinates  . .  5 

3.  Orientation  Angles  ^*1  Howabl*  Seer lom.  .  to 

*».  Variation  of  with  Fuel  3elgfat  Far».^eter»  a 

Cylindrical  Tank'' .  .  70 

5.  Variatioo  of  B#4^  with  fuel  Height  J’arftjester,  e 

Cylindrical  Tank .  71 

6.  Variation  of  Cylindrical  Frequency  Parameter 

with  fuel  Height  Parameter .  ?2 

7.  Variation  of  C,+j  with  Fuel  Height  Parameter,  t 

Spherical  Tank .  7< 

8.  Variation  tf  with  Fuel  Haight  Parameter,  e 

Spherical  Tans:  . . . .  .  74 

9.  Variation  of  Spherical  Jva^jttO'y  Paramter  with 

Fuel  Belgjt  ParaaeWr  .  75 

10.  Real  and  Equivalent  Retes^'.ar  Tank*  .  .  ...  ”?7 

11.  Reavltant  acceleration  tx»  Cylindrical  Tank .  8i 

*d.  Hcriioatai  Cylindrical  Tank  and  Spoerlcal  Tank .  Si 


A;" 

h 

Aj 

Cfc 

c- 

4 

c5 

E 

er 

F 

Fa 

**V 


SYMBOLS 

static  nerDcy:)fj3ic  'eras.  Se»  equation  (19S '• 
given  raoUe  of  deration  In  degree  of  freedom  *. 
component?  of  *  «,■  tr«*  y  coord  ii*’* te  syrter 

inertia  terrs.  See 

aerodvnaaic  stiffness  tersa.  Joe  Equation  (196). 

the  dyosslcaily  bt.lar.cing  rotation  rate  with  respect  o  q^ 
of  ta«  vehicle  relative  to  the  vehicle  axes. 

coapor.ent3  of  £j  in  the  y  coordinate  «v»ten. 

aercvlynaaic  daxplrg  terra.  See  Equation  (197). 

t.^e  UynaalcaLly  balancing  translation  rate  with  reaped  to 
qw  of  the  vehicle  relative  to  the  vehicle  axes. 

coeponer.tc  of  in  the  y  coordinate  systea. 
perautatl*n  avabol .  See  text  preceding  Equation  ^’39). 

nuaber  of  thrust  vectorirg  notzlec  (or  "engines"), 
components  in  the  j  *y*t*<a  n*“  the  3sj  vectors, 
the  *ua  of  the  external  forces  exerted  on  the  vehicle, 
the  external  force  on  the  h-th  particle  of  the  1-tn  section. 

Internal  force  charted  on  the  h-th  particle  of  the  i-ii- 
section  by  the  J-th  particle  of  the  Jc-th  section. 

the  nagnltude  of  FjfckJ* 
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HU; 

Hjk 

h; 
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toe  ivfi  of  thf>  nwats  about  the  origin  -f  the  /ehiclc  axe.* 

due  to  the  extern*!  forces . 

the  force  per  unit  ca:»  due  to  gravt*-*. 


the  coefficient  oi  ‘structural"  damping  associated  with  th* 
J-th  itgree  o:  frc.?dco. 

products  of  inertia  of  si  uc U. ~c  and  fuei  vehicle  axes. 

the  moments  and  the  negatives  of  the  prcduc.  •  jf  inertia  of 
section  1  about  Its  cvn  axes.  See  Equation  ''}8). 

inertia  coupling  tern*.  See  Equation  (173 )- 

modai  unbalances.  See  Equation  (150). 

subscript  -sed  to  donate  a  particle  of  t  section. 

the  partial  linear  velocity  with  respect  to  q^  of  the  center 
of  saas  of  section  1  relative  to  the  vehicle  axer  -  values 
obtained  after  dynaaic  balancing. 

ccsrporcnts  of  hj*  in  the  y  coordinate  system. 

acnenr-i  and  negatives  of  products  of  inertia  of  structure  and 
fuel  about  vehicle  axes.  See  Equation  (IU5). 


[  subscript  used  to  denote  a  section  of  the  vehicle. 

three  unit  vector*  pointing  respectively  in  the  directions 
of  the  three  vehicle  axes  yr.  See  Sec.  2  and  Pig.  2. 


,}  (,  three  unit  rectors  pointing  respectively  ic  the  directions  cf 

“  the  three  axes  *\)(  of  section  1.  See  Sec.  2  and  3  and  Pig.  2, 


the  partial  linear  velocity  vixh  respect  to  qk  of  the  center 
of  mss  of  section  1  relative  to  the  vehicle  axes  -  arbitrary 
values  given  prior  to  dynaaic  balancing. 


components  of  in  the  y  coord^ot#  systea. 


suffix  used  to  denote  e  degree  of  freedom, 
ccepooents  of  the  ctiffhes*  tensor. 


vii 


See  Zquatlou  (90). 


suffix  J  ?d  to  denote  i  degree  of  frccaoc 
aodaL  irertia  tcrrs.  Sec  Station  vl5I  ' 
suffix  u^ed  to  denote  j  j-gree  of  frc-dcc. 
corponcnts  of  the  terse-  See  Eq-ation  (32') 

the  bending  Kvr;  tt  a  specified  location, 
components  of  M  in  the  >  coordinate  rystee. 
total  aaci  of  vehicle  and  fuel  at  any  inrtant. 
cas*  of  section  1. 

nasi  of  the  h-th  particle  of  section  1. 
the  nuaber  of  sections  and  tanks. 

the  generalized  forces  associated  vita  inert!*  forces.  See 
Equation  (1*6). 

eodal  inerti*  terns.  See  Equation  (172;. 

nuaber  of  elastic  degrees  of  freedoe. 

a  unit  rector  located  at  a  certain  point  on  the  surface, 
perpendicular  to  the  surfece  at  that  point,  and  pointing 
ou'vard. 

ccepcawts  of  n  in  the*VX  coordinate  syatesa. 

the  generalized  forces  associated  *lth  conaervail.-  internal 
forces. 

position  vector  locating  the  origin  of  the*u£  coordinate 
systea  vith  respect  to  the  y  coordinate  systes.  Sec  fig.  2. 

components  of  o,  in  the  y  coordinate  ayates  ond  coordinate* 
of  the  fenter  or  ease  of  section  1. 

the  masher  of  particles  in  the  i-th  section  or  tank  of  fuel. 
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the  generalized  forces  asso-  luted  wit:,  d»  .si,'«tive  internal 
forces. 

»odal  accents  and  negatives  of  product*  of  Inertia  o!  Vx 
vehicle.  See  Equatic..  (S3). 

nodal  ncaent*  end  negatives  of  products  of  irertia  of  j*  t*  w. 
i.  See  Equation  (179). 

the  pspe*vl  '.red  forces  M^c-lcted  vito  ey  •inii  force*.  See 
Equation  (**7). 

generalized  coordinate  associates  with  t hr  w-*h  degree  of 
freedom 


suffix  denoting  the  r-th  coordinate  axis  in  eltnor  the  y  o- 
thc  "wf  system. 


the  surface  of  the  i-th  section. 


< 

suffix  denoting  the  »-th  coordinate  axis  in  either  the  y  or 
the  *v£  cysten. 


the  kinetic  energy  of  tha  vehicle  and  fuel. 


the  nagntldue  of  the  thrust  force  at  the  i-th  nozzle. 


ti»e.  Also  used  aooetijses  as  a  suffix  In  the  sane  sense  as 
r  or  s. 

potential  energy  due  to  elastic  deferaation. 
energy  dissipated  thru  Jarping- 

linear  velocity  of  the  venlcie  at  the  origin  of  the  vehicle 

axes. 

cesponenti  of  V  is  the  y  coordinate  aystea. 
velocity  of  th>  fc-tfc  particle  of  the  1-th  jectlon. 
the  work  done  b>  the  external  fcrces. 

the  "piston  spved*  (or  dnv^vath)  at  a  point  on  the  surface. 

position  vector  of  the  vehicle  In  relation  to  a  apace-fixed 
fraae  of  reference. 


position  vector  of  tha  h-th  particle  of  tha  1-th  section  in 
relation  to  the  vehicle  axes. 


cc?pcn-nt*  of  y^j  --  y  coordinates  of  the  b-th  particle  of 
the  ith  set.  tiers. 


U 


position  sector  of  t*^  c«**jter  of  ease  of  the  vei  icle. 


£ 


component*  of  yc  la  the  y  coordinate  system. 


4»- 


the  partial  S3gu.Ur  velocity  vith  respect  to  of  the  ni£ 
coordinates  relative  +'  the  y  system  -  values  obtained 
after  dynamic  calaccl:i£. 


components  of  OC  <n  the  coordinate  tyate«. 


the  partial  angular  velocity  vlth  respect  to  of  the  -of 
JV  coordinates  relative  to  the  y  syste-  art**'  — v  values 

given  prior  to  dynes!:  balancing- 
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component*  of  ,3^  in  the  "vj(  coord li r  < ' 


products  of  Inertia  of  sect*vo  1  referrt.  the  sectlooal 

axes- 


the  distance  from  particle  kj  to  particle  lh.  dee 
Squat ion  (77)* 


inertia  coupling  tana.  Sea  Equation  (171)* 

the  KronecJcer  delta 
§n  •  1  vhen  r  -  s. 
gri  •  0  when  r  /  a. 

the  logarithmic  decrement  associated  vlth  the  J-th  degree  of 
frtedoeu 
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*05 


Sc*  See .  3 . 


5. 


the  angle  or  rot»*lcn  o  1  3gi  aM  3’ l 


<r- 


ioe-tia  coupling  terns.  See  Equator  'l?o). 
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n- 
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o-ror 
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inertia  coupling  teres.  See  r*.  aticn  (l69)> 

nodal  product*  of  inertia  of  section  1.  Sei*  Kuation  (1^*2). 

asgle  of  swivel  of  nozzle  (or  the  J  ,  veruir)  “•>out  an  axis 
(\~)  perpendicular  to  3*  «nd  sakLig  -in  engl-^  .  with  3->* 

See  Sec.  9-  3 

Inertia  coupling  tern*.  See  Equation  (167). 

aerodynamic  *odsl  Urn.  See  Equation  (191 ). 

ratio  of  circua ftjenee  to  dicoeter  of  a  circle, 
ti-'*  a tso spheric  density. 

the  partial  linear  velocity  with  respect  to  q*  of  particle 
n  relatio*  to  section  i. 

cccponents  of  0^,^  in  the -\j^coordisete  system. 

position  vector  of  the  h-th  particle  of  the  i-th  section 
relative  to  the  origin  of  ,je -^coordinate  system 


cccponenta  of  in  tne“V(sy«tea. 
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angle  of  rotation  of  the  axis  and  plane  of  swivel  about  the 

y1  cxis  (3,).  See  Sec.  9. 

angular  velocity  of  tha  vehicle  axes. 

cooponcnts  of  A  in  the  y  coordinate  system. 


U)l  vibration  frequency  associated  *ith  the  3-th  degree  of 

freedom  See  Equation  (91 )• 


inertia  "symbols*.  See  Equation  (59)' 
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2.  IMPODUCTSP* 


T.*J «  report  include*  the  derivation  of  the  eqp&t'  .as  to  be  used  in  toe 
SUveti'isl  lo eda  Profiz<«  (SLP).  This  program  is  to  oe  used  >a  conjunct!.,*} 
v-ito  tr—  *lc  8Jjc-D*gr*o-or-7r*edoa  Tligfct  Peth  toefuter  Progra*  (SET),  a* 
a  —as  to  determine  the  reticle  structural  In***  sod  response  due  t<»  sero- 
oyrt«ic  leaf.-*.  loads  das  to  ccatrcl  surface  deflections,  sod  icjvizvnancttl 
ii.turbsrces  (l.«»*  riai  profiles  sni  ecs£is#«cus  dlsc^U:  t— butene*  profiles). 
The  proeraa  permits  the  taclnsjrn  of  up  to  i*i  elastic  degree*  ,-f  freyloe  ad 
40  fuel  slosh  aodee.  The  elastic  dsgress  «f  freedom  ur*.  •tMts-ry,  sad  the 
user  aay  Incorporate  say  mate?  of  sodee  such  ss  t*ody  bendiot,  viua  bending, 
wing  torsion,  etc. .  so  ss  to  total  17.  The  fuel  slosh  nodes  ^'orporete  2 
longitudinal  and  2  lateral  wiu  on  each  teas  sod  thr  pvofrea  allows  one  to 
Include  up  to  10  teaks.  It  is  recognlted  that  the  gross  vehicle  wtioo 
(lerc*  aotime)  Influences  t*m  easll  Motions  (elastic  deforest! ccc  sad  fuel 
sloshing)  of  the  rshicle£  text  It  is  sssunsd  that  these  sasller  aotloae  here 
%  negligible  effect  on  the  larger  Motions  of  the  vehicle.  Olher  basic 
atsuapUoo*  used  la  this  analysis  are: 

1.  Qndaaped  free  vibration  aodee  are  used  to  specify  the  elastic 
deforastions  and  fu«l  slosh. 

2.  Thera  is  no  elastic  cr  darpdeg  coupling  between  the  degrees  of 
freedom, 

3.  The  aerodyasnle  forces  esa  be  obtained  by  ftcvtcnlan  flow  theory. 

k.  Tho  fuel  surface  (excej  'or  the  sloeMng)  is  cooeiderad  to  be 
perpendicular  to  the  rt  .tart  aeeelaraU.cn  at  ths  center  of  the 
tank. 

5.  The  fuel  sloch  nodes  of  c  tank  that  is  sot  vertical  or  larlsootal 
tan  be  rvpraewted  by  those  of  sons  hypothetical  task  that  Is 
vertical  or  heritor**!.  Ir  addition,  longitudinal  fuel  clashing 
In  s  horlsoatsl  eriindrlcal  tank  is  reprae  anted  by  aa  analogy  to 
a  rectangular  teak. 

6.  its  effect  of  a  rocket  eogln*  esa  be  reprsr acted  by.  a  \  irust  teeter* 
idach  is  a  si^Ufitetion  that  ssi.nea  the  center  of  ossa  flenr 
through  the  aossl*  to  be  exactly  allusd  with  the  geometric  axis 

of  t be  noetic. 

The  coaglexltles  iahsraat  la  this  type  of  pcohlsa  ax*  so  great  that 
certain  conventions  of  the  tensor  notetim  are  iaecrporated  in  the  xubsoq^rt 
development  in  order  to  shorten  the  writing  of  the  sections.  These  opera¬ 
tions  ere  explicitly  explained  ae  they  are  Introdnced.  The  soalyeie  of  the 
structure!  loads  is  loyally  developed  in  the  following  sequence: 
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Tv'“  v-hiclf*  kiueuaties  arc  derived. 


2.  Tb-i  force  and  moment  relations  ore  found  using  Vtvtonlan  mechanics. 

3  '’’he  equations  of  notion  of  tbc  c2aatic  deforuatioc  are  derived  (work 

*»r:d  -nergy  concept#  are  used  J.o  check  the  basic  iormui&tione  of  Items 

c.  tu.'i  3.). 

U.  ihe  &ain  equations  to  be  »••«*  to  determine  the  elastic  deformations 
are  put  into  tertn  suitable  for  computation. 

5.  The  aerodynamic  forces  (using  Newtonian  flov  theory)  *i-e  found. 

6.  The  analysis  of  the  fuel  >lo«h  problem  i«  included. 

7.  The  thrust  forces  are  introduced. 

8.  The  accalerationa  at  all  locations  are  found. 

9.  The  shear  forces  and  bending  moments  ere  calculated. 

The  generalised  forces  to  be  used  In  the  program  are  Inertia  forces  Njf  exter¬ 
nal  forces  Qj,  conservative  internal  forces  Oj  and.  dissipative  internal 
forces  Pj  .  These  forces  are  represented  by  Equations  (46)  -  (4?). 

To  clarify  to  some  degree  the  subsequent  analysis,  the  representation  of 
t  *  coordinate  system  ij  presented  in  Figure  1.  The  origin  of  the  orthogonal 
reference  frame  (  >4’’  •4'  )  is  represented  by  an  arbltrar;  point  that  would 

be  fixed  in  the  vehicle  if  it  remained  rigid  during  the  notion  along  its 
flight  path.  In  conjunction^ with  this  frame  of  reference,  ere  located  relative 
coordinate  systems  ( l/  ,~L>‘  ,  X)  )  positioned  at  various  points  on  the  body 
to  define  the  elastic  deformations  and  fuel  slosh  motions.  An  absolute  refe¬ 
rence  frame  ( ,  'A  ,  X*  )  is  shown  for  gensmllty  with  *  the  position  vector 
connecting  the  origins  of  the  reference  frames.  As  a  physical  insight  into 
the  relative  relations  of  these  coordinate  systems,  cons,  dor  the  rare  when 
the  vehicle  center  of  gravity  (C.O. )  is  the  wtlgln  of  the  vV  .<-*■  ,  M-  triad; 

then,  the  velocity  of  this  point  it  represented  by  >Tc4f  .  6 
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Figuu  1.  Coordinate  Systems 


■  KIiatATICb  Or  TBS  VEHICLE  AJfD  BASIC  ttSSttHPTxOSS 


la  JJV.lyilCLC  the  Button  of  \  vehicle  In  flight,  It  Ii  r  avenlent  (peri-apa 
occessar-<  u>  u..r«  In  tern*  of  the  fo Having  type*  cf  Botlc”:  (,l)  the  notion 
of  the  vehic*v  as  *«  whole,  vhlch  charactex ices  its  ■’flight"  and  is  referred  to 
here  t UK  action  cf  the  vehicle,  (2)  large  •otico*  of  vervain  parts, 

•wch  iuj  ccatrol  surfaces,  relative  to  tb*  reut  of  the  vehicle  -wad  large  dle- 
rloccse.it".  of  t>-c  fuel  -c  the  tan*t>,  an!  vy;  s nax-  <lwui  t**i«r*eti  one 
r.cl  sloshing.  Types  of  notion  (l)  and  (2*  are  deteradce-1  In  the  basic  Six- 
Lesre"*cf  Freedom  Flight  Path  Study  Generali  ted  Computer  Progra*  ( i>LY j  and 
the  Vehicle  Fbyslca'  Characteristic#  Subprogram  (VPCS*.  Tyem  (3'  Is  to  be 
determined  In  the  Structural  Loads  ?rogrc»  (SLP),  vhi tn,  as  ltr  nvse  inllcates. 
is  also  to  acteneine  the  •  U-uctural  loads. 

It  is  assumed  here  that  the  relatively  snail  elastic  doforsaUnae  and  fuel 
sloshing  Motions  have  a  negligible  effect  on  the  other  (large)  notions  of  the 
vehicle.  (Some  considerations  associated  vith  this  assumption  are  Investigated 
In  following  paragraphs. )  It  ie  not  assumed  that  the  large  notices  of  the 
vehicle  have  a  negligible  effect  on  the  snail  notions.  Consequently,  the  large 
notices  *  types  (l)  and  (2)  -  will  be  employed  ae  part  of  the  input  to  the 
Structural  Loads  Program. 

An  orthogonal  right-handed  triad  or  unit  vector*  J,  ,  J4  ,  J *  that  would 
be  fixed  in  the  vehicle  If  it  vex*  perfectly  rigid  1*  introduced  to  provide  a 
franc  of  refererce  (a)  to  represent  the  gross  notion  of  the  vehicle  and  (b) 
to  facilitate  the  description  of  the  Ojbe^  notice*  of  the  vehicle  -  types  (2) 

*rtu  (3)-  Rectnnguiax^coqjdinatcs  'i  '"i  ,xt  ***  associated  respectively  with 
the  unit  vectors  J,  ,J«,  J*  au  shown  In  Figure  2.  These  coordinates  sre__  t 
seen  to  be  the  conponents  of  the  position  vector  ,  which  equals  J, ♦ 
The  axe/?  of  these  coordinates  are  called  'vehicle''  axes.  * 

Additional  orthogonal  right-handed  triads  of  unit  vector*  J  .  «0«  •  J* 
that  would  be  fixed  in  the  various  parts  of  the  vehicle  end  in  the  fuel  in 
the  various  tanks  if  they  were  rigid  are  introduced  as  frsnet  of  reference 
(a)  to  represent  the  notions  of  the  perts  end  the  displacements  of  tbv  fuel 
relative  to  the  vehicle  and  (b)  to  facilitate  w*  description  of  the  elastic  (> 
deforest! cc*  and  the  sloshing  of  the  fuel.  Rectar^U*x ^coordinates  V  .  V,  V 
are- associated  respectively  with  the  unit  vectors  J,‘ *  J*  *  J*  and  these 
coordinates  arejthe  components  of  the  vector  ,  vhlch  is  the  position 

vector  in  the  J '  coordinate  system.  The  exes  of  these  coordinate*  are 

called  ’■section"  axes.  The  vector  O  locates  the  origin  of  the \  coordi¬ 
nate  systen  with  respect  to  the  Jr  systen.  Ccnrequently,  then,  •  S'  . 

The  gross  notion  of  thejrehlcle  is  that  of  the  JrCv*  __trl*d, 

viuch  has  a  linear  velocity  v  at  its  origin  and  an  angular  velocity  SI  . 

These  velocities  are  functions  of  the  tine  t»  ,  and,  together  vith  their 

derivatives,  completely  describe  the  gross  notion  of  tbs  vehicle. 

Generalized  coordinates  are  employed  to  specify  the  cugflguretion  of  the 
vehicle  and  fuel  relative  tc  this  frsne  of  reference  -  the  Jr  triad.  In  so 


Figure  2.  Vehicle  and  Section  Coordinate* 


doing,  However,  *  distinction  Is  cede  between  the  *arge  not of  type  (2) 
and  the  tneil  notions  -  type  (3)  •  tnaaxuth  e><  the  iarg-  uot'.or*  are  foreknown 
in  the  5II*  they  can  be  specified  by  neons  of  a  jingl*  .;cr*raiized  coordinate, 
for  Milch  tae  syebol  cf  is  chosen  here,  and  vfaich  is  tw  equated  to  the 
tire  f  Tn..  use  cf  Cf*  for  this  purpose  rather  than  r  ,  even  t house  the 
tvj  --  r' lyrically  equivalent,  serve?  *c  <*ifitiagu*  »n  th%  notions  ot  tvpe  (2‘ 

tv*  motion,  type  (i),  the  gross  notion  beiiv  r.-prccested  as  a 

function  of  t  but  not  as  a  function  of  <j° . 

Uttcasipad  free  vi oration  sodas  arc  used  as  degrees  of  fruedos  for  sjecify- 
elastic  deforest tors  fuel  sloshing,  tactions  of  type  (3),  and  ar» 
referred  to  aa  elastic  degrees  of  free den.  (Gooj  results  can  reasonably  be 
expected  if  a  sufficient  nuscer  of  the  lover  frs q— acy  sre  us*d. )  These 

stall  notions  ere  specified  by  the  generalised  coordinates  "  ,  o\  ..  o*  (  q 
being  the  nuc.be r  of  elastic  degrees  of  freedoa).  *  '  '  ' 

At  this  point,  it  is  convenient  to  ■»dopt  the  range  and  floatation  conven¬ 
tions  of  the  tensor  analysis  as  follows : 

(1)  Range  Convent  ion  -  A  coordinate  suffix  that  occurs  Just  once  la 
a  tern  is  understood  to  represent  ell  the  integral  values  appro¬ 
priate  ’to  its  range. 

(2)  Suaaatlon  Convent ion  -  A  coordinate  suffix  that  occurs  Just  twice 
In  a  tore  ioplies  suaaatlon  with  reapect  to  t  vst  suffix  over  its 

range. 

Yhev  ocssvsitions  enable  us  to  write 


5 » jr  <f  (=j.v  vM*‘  a) 

i;  =  J'  ~j'X=  J.V*. JjV*. J iV*)  (?) 

5  =  Jr»' \(f-j_,e'+  oJJ  (3) 

J*  &r  t- J, e  r  4  Jg.  €.  t)  (* ) 


H»«  &r  in  equation  (3/  «re  the  coordinates  in  the  Jr  coordinate  systea  of 
the  origin  cf  the  J'r  systee.  The  e'r  in  (fc)  are  the  cocpooenta  in  the 
ays  tea  of  tbs  J’r  vectors;  for  any  particular  choice  cf  r  and  5  , 

(?r  =*  J\  -J  $  ,  which  is  the  cosine  of  the  angle  between  Ji  and  J s  , 

Equations  (l),  (2),  and  (3)  illustrate  the  use  of  the  suzeatlon  convention; 
equation  (U)  illustrates  both  conventions. 


Unless  otherwise  noted,  the  range  of  the  suffix  of  a  unit  vector  (  Jr 
or  Jp  )  or  of  a  rectangular  coordinate  (  yr ,  *v'r.  or  9r  )  is  1,  2,  3. 

The  range  of  the  suffix  of  a  generalised  coordinate  (  <?c  )  will  be  understood 
to  be  !  y  Z .  •  •  •  n  •  Zero  (aa  is  )  is  specifically  and  deliberately  -?>• 
eluded  in  the  use  of  the  range  s*sd  suoaatloa  conventions.  Thus,  In  specifying 
basic  functional  relations,  vc  write  „ 
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vt?  find  by  substitution  fro*)  (!'  thr  (*i  that 


—  r  —  r  r 

j.u-.lrc  r  i.  v 

_  v  '  -  >  ' 

'  Jr  O  *  J^dr  'J 

(9) 
(10) 


jr  V 


6  *■  J  =  a  t' 


The  components  yr  of  the  position  vector  y  are  the  coor'4 Lea* es  of  a 
i^rticle^of  the  vehicle,  and  their  a*  tuc  pcrticle  moves  vlth  respect 

to  the  Jr  frame  of  reference  is  a  function  of  and  toe  O*  .  oince  the 
angular  velocity  of  the  Jr  triad  is  /l  ,  ana  since  the  »~r  /arc  functions  of 
the  tioe  t  only,  their  cerivatives  are 


~di^^  ^xJr 

It  is  clehr  fros  (10)  that  —  O  ;  therefore,  frtw  (1), 

^  U  -  j  jr  j'  -  T\  i  J .  L{  -  C\  t  U 

3f  <1  C  J  J  J 

’iff 


i  *  U 
J 


since  O 

/ 


£oc . 

it 


in; 


(is) 


(13) 


(11) 


and  is  unknown  anti.  Jet  trained  in  the  so*  lion  of  the 

equations  of  notion,  which  are  to  follow.  Man  log  use  o;  *  id  the  fact  that 

the  lixy-**r  velocity  of  the  origin  of  the  Jr  tnec  it  V  ,  ve  find  that  tne 
veloeitv  of  n  ^.article  of  the  vehicle  is 


<vi —  -  1  /  ^ 


j  f 


(15) 


The  acceleration  of  a  Dart i cl e  can  be  found  at  follo>* 


Likewise 


ft  (0) =  ^  *  If1  '  <5^V  +  '5^4"*  ? L  • 


U6) 

(17  > 

(18) 


The  acceleration  is  now  round  by  differentiation  of  (15)  and  substitution 
froo  (13),  (IT),  (18),  with  the  result 
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slO  -  T  dlL  +  sLL  O' 

,-jr  r  H '  dx 

=  J,  i  t  -f-  /I  r  -»  l 

-Sr^, 


'?0j 


A  *  Jr  l\  =  At  y  A  =<P. 


The  .A  -  ♦  hen,  are  the  coapcnents  of  the  angular  accelerate::  The  linear 
acceleration  M  the  origin  it 


~jrv r +n  *  v 

tk(vl+fiV-  A'/) 

ij3  (v3*  nv*-/iv) 


(21) 


The  coefficients  of  the  unit  vectors  arc  tne  cosponents  of  the  linear  accel¬ 
eration. 
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Ta  IV  application  of  *evton*s  sec  cod  law  of  aotJ^,  to  the  /eruci.-..,  i- 
is  oc *-  have  a  •seans  of  identl fyirg  the  particle*.  But  the  -efi’-ie 
i-  Into  varies*  part*  (or  sections ^  and  various  fuci  taoxs,  vhi;c 

*l*o  c*~*  *f'  oe  identified.  Because  of  It*  shifty  and  siooby  net-r*#  toe 
-x-el  Coiv^i  be  regarded  a*  part  or  the  tank  that  contains  it.  The  tan*  it¬ 
self  i*  treated  a*  oce  or  nore  structural  st.«  cna-  Suhner* pt«  ere  <*“i  .yed 
to  i(Unbif>  j&stes  or  «ass  particle*  and  their  rectangular  coordinate*,  a 
■'ir-^ie  t  bs^rl z*  or  *h*  r«rst  t»»o  \»e  vectio*  or  the 

fuel  contained  io  a  certain  tank,  and  the  second  subscript  d* noting  the 
particle  of  the  section  cr  fuel.  The  absence  of  *urn  £jt*c  ij.la  denotes  a 
quantity  pertaining  to  the  entire  vehicle. 

Thus  the  a**^  of  the  h-th  particle  of  the  i-tn  -<rctioa  zc  /r.(  ^  and  it* 
coordinate*  are  U  th  and  v\h  •  The  oass  of  the  i-tb  section  is  m,  and  its. 

** coordinates"  a&  &{  .  The  oast  of  th*  entire  vehicle  i*  m  l  without  a  sub¬ 
script).  Let  f\  be  the  nunb«.r  of  particles  in  the  1-th  section  or  tank  of 
fuel  and  If  be  the  nuesfcer  of  sections  and  tanks. 


then 


and 


/r 

rr  - 


rri 


(22) 

(23) 


Let  the  Jr  triad  of  the  i-tb  section  or  tank  of  fuel  be  designated  as  the  3*. 
triad  and  let  its  origin  be  at  the  center  of  taass  of  the  1-th  system  of 
particles.  Then  the  a'txrc  the  coordinates  of  the  center  of  sas*  of  tbc  i-tfc 
section,  and 


A 

21  «»*  t\s -  ° 

/W  {2k) 

Also,  let  he  the  coordinates  of  the  center  of  *a*s  ui  *,*,e  v  hide.  Then, 
vlth  the  aid  of  (10),  (22),  and  {?'■*)  it  is  found  that 


Pi 


3c 


A- 

•!—  2__  r".t  3  I* 

l»/  hlf 


-- £  t  ».  (< * <, 

-  ZT «; 

»*/  (23, 
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w3ition  vector  f  the  center  of  a&cc 


-  .Ui  tl*  '.id  of  Newton's  1«.-  of  —  ‘.zr.,  1.  ~*r  _h<,  vn  tha*  tit. 

k  *  •?  forces  ex.rtcd  on  all  the  particles  oi  t**e  v*f-icle  cqua  s  t*i 
su~  of  tb*  ' 'terra*  forces  only  tot  this  r .  lesigyiaied  by  f,  c.i 

"s;o"i  1a*  -f  -notion  cod  'If), 


»»/  n-4  tt  erf 

^  It  <-,|#  ’ '(“'.y-v f  fp" '  y.* 
,i(ri 

4%f] 


(27) 


Likewise,  the  sue  of  the  ootsents  about  the  oriji.-  of  the  Jr  triad  due  to 
all  the  force#  la  equal  to  the  nua  of  the  sooenta  due  to  the  external  forces 
only.  Let  this  oe  designated  by  5,  then,  by  Newton's  second  law  of  rcotion 
and  (19), 
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Incsmuch  ns  the  Jr  triad  provides  a  frame  o'  reference  ■"‘present  the 
£ro_ps  ration  of  the  vehicle,  its  linear  and  angular  velocities  via  acceleriticn 
V,/\,  c”/dt,  and  u.//dt  are  those  of  the  vehicle  as  a  whole,  it  w^old  be 
^trictiv  proper  to  require  these  velocities  and  accelerations  to  satisfy  ( 27 ; 
and  (2  o', ;  but  H  has  been  assumed  that  the  elastic  defo'sations  and  fuel 
sloshing,  rations  of  type  \S)j  have  a  negligible  tffec.1  on  the  large 
of  t.h«  v"h '-itt,  therefore,  Jgi,  dV/dt,  and  dtydt  are  -egarded  as  not  being 
fuiu  * ‘.'.-"to  of  the  <rene;-alized  c.-.ordinatca  /  or  their  derivatives  )*and  *’* 

The  fe  •*  ti-‘,  P  and  6  may  b„  significantly  affected  by  the  elastic  deforma¬ 
tion  oi  aerodynamic  surfaces  is  arbitrarily  disregarded  here,  and  the  portion; 
of  (27)  and  (28)  involving  y  “  ana  y  are  simply  ignored  in  t!:<:  process  cf 
determining  the  gross  motion  of  the  vehicle.  This  leaves  ",  for  the  determi¬ 
nation  of  the  gross  motion, 


/:zt  vr  "’••{df 


iih  t 

8  <1  J  7  >  7  -» 


129) 


G 


v- 

rr  £t  " 


d?  *•  A  >  (f--  >  <ju)  ♦  {$-  xy.< 


•■J  A  x  ayu.  f 


(30) 


with  none  of  the; e  terms  being  regarded  as  functions  of  the<^‘ 

*  ' 

Equation  (2b  j  contains  the  summations  ST  £_  m  .  K  h 

,.i  >\.i  •"  >'‘n 

v  r  . 

and  21  £  /nif(  jkfr-  .  If  the  terms  of  (29)  are  not  to  be  functions  of 

the  y  %  the  partied  derivatives  of  these  summations  with  respect  to  the  y* 
should  be  equal  to  zero.  Furthermore,  if  these(saoe  partial  der.  .tives  are 
exiual  to  sero,  the  terms  of  (27 'I  involving  cue  y*  and  the  y*will  vanish, 
because  they  contain  tuese  partial  derivatives  as  factors.  In  fact.  It  is 
sufficient  for  thin  purpose  for 


to  be  zero,  because 


*•/ 


f; 

V 

4-  m, 


#7* 


if 


J^jr,  m  IZ  rr  IZ  Tl  -  O 


v  -x. 
V 


r,  £  -"..fSjt-a 


Reference  to  Equation  (26)  sheds  a  little  more  light  on  this  problem: 


£  If*  =  ^  =  x”  Ifr- 


(31) 
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Frca  tJ.e  physic* 1  viewpoint  it  is  clear  tnat  {2  will  not  rhaage  uuca  ns  «i 
result  of  elastic  de formation,  but  that  fact  does  not  preclude  the  po  slbinty 
of  its  -iMwiglng  rapidly,  therefore,  it  is  rash  to  recuse  arbitrarily  tiat  (31) 
will  equal  zero.  Satfcer.  it  is  de*i  »-able  to  lapose  its  bring  zero  as  *  ccrii 
tion  w‘f  satisfied  by  the  elastic  degrees  of  freedom. 

Ta-  :oapcncnts  of  *he  inertia  tensor  (or  eleaents  of  tr.  inertia  -satrix) 
of  U. e  .*  \e  are  given  by  the  formula  (Reference  7,  Equation  2-12) 


(32) 


It  is  xncvn  that  vhes  the  J-th  degree  ox  ireedee j  is  a  for-  of  soticc  ia 

which  the  rt-hic)**  ’ccv,e*  ***  •  rijii  body,  laving  ®  trvnslation  rsue  C ^  and  a 
rotation  rate  t>j  relative  tc  tbe  j  r  triad)  and  tht  ts-th  degree  of  freedoa  Is 
a  noraal  free-free  »le  of  v.bratioo.  Vhen  such  is  the  case, 


and  froc  (32) 

(c>  *  y ' 

8  t  ir  "off  4 Jj  •  £  ^ 

=  0 


(33) 


(3M 


flow  Tj  and  ,  being  any  translation  and  rotation  rates  of  the  vehicle,  aro 
arbitrary;  therefore. 


tt  mu-^r:0 

m4  s  , 


(35) 


(36) 


vhsn  the  k-th  degree  of  freedom  is  a  norool  free-free  lode  cf  vibration. 


I 

| 
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Th«»,  if  the  elastic  dcsre*»3  of  freedoa  satisfy  the  cor^L. ;1  'n  ihat  the/ 
are  .  ,rcdi  ,'t:  frc.  nrdcc  of  vibration.  Equation  { 35)  is  oatis-»rd,  the  i-n*. 
la  (2T)  involving  <t‘and  n K  ’  z^‘  rh  (c’ia»oavins  the  differe*  ,e  between  [2 1, 
and  (2>J),  and  the  terns  on  the  right  side  of  (2$)  are  mt  furctions  of  tne 
<a*.  However,  the  use  of  nor=al  free-free  nodes  .*.co  ipu-shes  oore  than  this, 
equities  ( 3w.  directly  el-1  f*1 natea  one  tern  of  (2d),  differtntiaM~-  ,f 
(30  ieaoj>  wo  the  dial  nation  of  other  ter**,  as  f ol  I  *  -• 


=  o 


m 

because  interchanging  the  superscripts  k  and  1  does  cot  affect  the  left  side 
of  (3?)  whereas  it  reverses  the  sign  of  the  right  side,  and  only  cero  equals 
ita  opposite.  Purtheraore,  the  superscript  L  could  be  replaced  by  C  in  the  tvo 
*  (justices  above;  thus,  (28)  is  reduced  to 


G-trh  +A 

nil, 

1  1  (39) 

which  could  be  used  for  whatever  value  it  sight  have  in  solving  for  the  o*,  it 
being  recognised  that  the  teres  of  this  equation  are  functions  of  the  4*,  in 
contrast  to  the  use  of  equations  (29)  and  (30}«  1 

The  results  thus  accomplished  by  the  use  of  normal  fr*»e-frw*  andes  of 
vibration  can  also  be  brought  about  b;-  a  "dyaaadc  balancing'  of  each  degree 
of  freedom,  individually.  In  order  to  do  this,  let 

-  a  ~tb  x  U  +  C 
TU*>  J*h  * 
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(40) 


the  Hc  being  given  (not  necessarily  free-free)  modes  of  vibration,  ana  the 
and  6 „  being  as  defined  in  connection  with  (33)  except  that,  lactead  of 
being  arbitrary,  thev  are  now  unknowns  to  be  determined  in  such  a  way  that 
(HS)  and  (36)  will  be  satisfied.  Substitution  from  (bO)  into  (35)  results  in 


V  •> 


(*...  +  K  *  u +  e.) 


and  cubstltution  into  (36)  results  In 

w  o 

H £7  "’ll  9,h '  (?«•>.  +  k« '  iU +  c«) 

=  tr  h  Ju.  *  (8„„.  ’K'S-i)*  «■£  1 


Mow  let  us  eli adnata  c,  by  forming  the  vector  product  of  Jcwlth  (bl)  and 
subtracting  it  from  (b2).  This  results  In  " 

IS™..  jik  *  (a«*+V&..) 

-  &  *££  ni„  *.<»  - -  o,  W) 

which  can  be  solved  for  the  .  Once  the  B  ,  are  otthlned,  (V* can  be  used 
to  obtain  the  c*.  When  the  br  oado*  ere  obtained  in  this  *..voer,  the  use 
of  (bO)  results  In  <3^-  that  satisfy  (35)  cad  (36).  These  asy  ha  called 

"dynamically  balanced"  nodes.  They  have  the  practical  advantage  of  being 
much  sore  easily  obtained  thua  the  normal  free-free  nodes. 
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h.  EQUATIONS  OF  MOTION  Ftft  THE  ET/.3TIC  IEFCRMATr~«<S 

I;.  the  preceding  section,  the  Influence  of  the  Intern*!  forces  and  the 
dtctr flutter,  of  the  aerodynamic  pressures  over  the  surfc.ee  of  the  vehicle  were 
ucllucVitcij'  d  •.  ..regarded .  In  this  section.  It  will  he  necessary  to  give  thesi 
full  cun * 1 vexation,  because  their  effect  on  the  elastic  deformations  cannot  be 
dls.-^rded  and  because  the  purpose  of  this  section  is  to  deduce  equations  of 
motion  tr«  tKC  determination  of  the  elastic  deformations  i;.  the  various  degrees 
of  frac-km.  For  the  sake  of  suitable  notation,  let  denote  the  external 
force  on  th-  h-th  particle  of  the  1-th  aeeti.*.,  and  let Tidi ••  represent  tine 
internal  force  exerted  on  the  h-th  particle  of  the  1-th  section  by  the  J-th 
particle  of  th.»  k-th  section,  fy  Newton's  second  lew  of  motion,  then,  the 
total  force  exerted  against  the  h-th  particle  of  the  1-th  section  i« 


m 


The  equations  of  motion  in  terms  of  generalised  ‘  >rces  are  obtained  from 
(44)  by  forming  the  scalar  product  of  ytt-d.  with  each  term  and  summing  ovtr 
h  and  1.  Thus 


N  Pi 


For  convenience,  the  generalised  forces  are  separated  into  four  types  and 
designated  as  follows: 

1.  Those  associated  vlth  inertia  forces  are 
to  Pi 


(45) 


2,  Those  associated  with  external  forces  are 


(46) 


Qi- 


(47) 


3.  Those  associated  with  conservative  Internal  forces  are  Oj. 

4.  Those  associated  vlth  dissipative  Internal  forces  are  P, . 


Since  tdie  f^£4^denote  the  Internal  forces,  ve  may  let 


i 
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a  no  sub-titi.t  ic/i  frox  tbefce  last  three  Into  (1*5'  leaas 


:«o) 


&) 


Substitution  -  r^=  (19)  into  (!■& )  results  In 

u  p. 


tf^TS  *-  a '  (fr'M+  JT 

+  2  n  -V  2(Qx^H^6k  (5o; 

+S?&  W  + 

If  we  cake  uae  of  (32),  (35^,  (36),  (38),  aad  acne  new  ayabola  In  «o  exeainat«oc 
of  the  individual  tenia  of  '50),  ve  obtain  a  siopler  and  cere  practical  ex¬ 
pression  for  H,,  aa  follows: 

V  ft.  ITT  ,T  *  * 


lp^-%  *S  •&£*»»“  Sjf- 0 

=  1 1  ru  an  •  f -)fri  •?*)-($*•  H)(^l 
=  f  l  a*  $£«■  -  cw  ^  ^Kj 


(5i> 


►‘1(52) 


D 

•  hfj 


t  r  ft.* |*r  ■  jf  -¥**$?  '|l|,ll,l-<-?1<',^i 


A. 
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f  fx.s  &-L  ■  =  foo-?! 

^iVi  0M  1 — — ■ ■=-'  (%) 

it  "•<  $*• -5  -  ^  -  a  •  i  |m.i 


=  O 


U,  •«•*$* -^.-s 

It**  ■  a 

li"u^-^<'-  h"  - 

Substitution  frcn  C?1 )  through  (60)  Into  (50)  results  in 

Nj  =  ~  CTQ?  P„, j  -rjoo.jj  +  a.  [ovC/i]  C\*- 

+£13^  $  *-  M  \*  cp .  (61) 

----  -  %  raw 

»a4  frca  (15)  en4  (IS)  that 

The  kinetic  energy  le  given  ky  the  veU  taova  fonula 

T=i£r*^  , 

Vhence,  vith  the  aid  of  (6?)  awl  (63), 

il.  ^ 


(5?) 

(56) 

(55) 

(t>o; 


(62) 


(63) 


(61) 
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-tz  wu  ^  , 

( £.}-  i  l  wu  i  + ^  i:  ( ^)] , 

tv  •  fi.,m-‘  ^  ■'St 

Subtraction  of  (67)  free-  (CC)  results  In  the  lA^rtn^  an  erpr_*/3ioa 

iW-ft'liSS*  ■£* ■  »■ 

by  the  defining  equation  (to).  The  use  of  this  expression  to  obtain  (6l) 
lead*  to  tie  interesting  discovery  that 

ft-nsp„i"  ia'ns 

Trf  1  t.j  *  P»vj  , 

i«s  ||W(s«fr?«.^-  #:*  *•**) 

-  accent*  and  negative*  of  products  of  inertia  of 
•t ructure  and  fuel  about  vehicle  axes. 

It  1*  also  int«r«*tli*  and  useful  to  obaerve  froo  (60 ),  (o2i,  (65), 
and  the  fact  tbit  *•  not  a  function  of  the  q*  that 


or  that 

where 


Thte  i«  eep*claily  u.tfMl  lr.  treating  tbe  Inertia  effect*  of  fuel  »loeb. 

Recalling  tb*  defibitlob  of  F.OU  '  w  loov  lee  nevtoa'i  third  lav  of 
sotion  that  0 


(67) 

(68) 

) 

(69) 

(TO) 

(71) 

(72) 


and  that  ^4^  is  parallel  to  j 


<73) 
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UAj.  =  o  (T- 

Substitution  free  (73)  into  (43)  results  in 

If  (4S)  in  idled  to  (75)  ud  the  sub  divided  by  2,  the  recult  ic 

oj  -  pi --tilil  r^t 

-  \u-w  -  a4: 

which  ic  the  distance  froa  particle^*  to  particltcil ;  and 

,rt  m^-j  •  a* ,  w 

positive  when  it  tend*  to  Increase  a^d  negative  when  It  tend* 

to  decrease  .  ' 

Then,  ilac*  TT^^ic  parallel  %  it  can  be  shown  that  (76)  ic 

equivalent  to 

Let  U  be  the  potential  energy  due  to  clastic  deformation,  sad  ■»•♦:  v 
be  the  energy  dissipated  through  dating*  Tbeee  both  represent  work  done 
in  overcoalag  internal  forces;  therefore, 

-  -  -kSxu!,  ,^V 

+  ^  (?**•"  ] 
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because  of  ^ 15 ) •  The  tera 


i\M.,  fl  >  ( h<-  /Vr  '  &  -  O  (81) 


o?c&”s~  «*■  tbersfore,  because  of  (76), 


Since  the  0*  ere  associated  with  eonaerratire  Internal  force*  ard  tb* 
P.  are  associated  with  dissipative  Internal  force*,  it  1*  clear  frea  tb* 
definitions  cf  U  and  V  and  froa  (82)  that 


Ifr" 

1  Oo-t-O)^ 

(84) 

“d  HL 
cit 

-  p.  *  f  . 

(81) 

Sow  U  Is  a 

function  of  and  tbe  Cj4  but  not  of  their  tl*e  derl^atlres; 

therefore, 

ft 

>■$? . 

Proa  (8h)  and  (86),  ve  see  that 

0, 

(87) 

Ibis  can  he  extreesly  helpful  In  the  cccputatlcn  of  tbe  0, . 

If  tbe  are  cclj  the  dissipative  lots  reel  fcrcec,  equation  (79) 

nay  be  directly  useful  for  the  calculation  of  the  P. .  In  us1e«  e«ptatlon  (47) 
to  cotQute  the  Q.,  It  la  not  necessary  to  Include  toe  force  due  to  grtvity 
because  such  a  force  can  is  represented  as  *n  v^c^‘  c*c* 
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wt v*n  1*  aatiafied. 

-s'  .s  introduce  certain  aswrcpticur  here  and  ..roceed  to  scoe  furt^r 
*?*+***•-+  c*  the  C.  and  P.  First,  let  us  aasuae  that  U  i»,  a  ainiaan  vcen 
the  cu  .-qu*l  tero;Jthen  J 


Oj^ij.O  «Ae-  «*•  ^  -  o  ;  (89) 

sad,  as  a  close  and  convenient  appnrriantlc©  (Reference  7,  l .  atlon 


<*■  Ki*c|K  , 

vhere  l<  4iV_ 
J'-jW  5o^». 


evaluated  for  the 


cf  »  0  . 


(90) 


Ina—ich  as  undamped  free  vibration  nodes  are  used  as  degrees  of  freed co 
for  specifying  elastic  deforaatlco*  and  fuel  sloshing,  there  is  a  frequency 
Uuj  associated  vith  the  J-th  degree-  of  freedcs  for  all  the  values  of  J.  For 
the  first  degree  of  freedom. 


or  K„-*  .  t?2) 


Likewise, 

W«. 

K« 


(93) 


and  so  forth  ft*'  all  the  degrees  of  frtedoa.  It  is  nov  further  aasuaed,  and 
this  rust  he  carefully  noted,  that  the  degree p  of  freedcei  vill  he  so  chosen 
that  there  vill  be  no  elastic  coupling,  that  is,  so  that 


K 


iK 


O  wkeh  i  7*  K. 


(90 
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oaring  tbt  ‘egrees  of  freedoa  in  this  'ashion  is  a  co'ao.i  practice  *n  tie 
analyein  of  i  i-jttc**  stability. )  A  general  expression  fo.  the  j;  t^ie 


= 


vh»*«»c».«  wtlcn  late  (90)  yields 


*95) 


0j  -  •  '  ,v>  * 

These  equation*  \“1*  tortKgn  (96)  kSt  based  on  the  cathaiatiCAA  relatione 
expressing  the  vib.  &to~y  motion  of  the  eysWa  in  c«u^  toe  given  degrees 

cf  freedcs  at  a  tlae.  A  further  pursuit  of  this  line  of  thc-uht,  linked  with 
the  association  of  a  coefficient  oi  r  structural*  d/cping  g  ,  w.  th  c-*»ch  degree 
of  f-eedcre  leads  to  a  oicplc  fora; la,  aialogous  to  (<£)»  fo't  /*< .  This  is 

!S7) 

The  determination  of  ».  free  the  logaritbaic  decrement  &'  is  sieple,  os 
follovs  J 

OT;.  a  _££j _  (98) 

«-Mir  *+($*)* 

—  vaen  6*^  is  ssoll. 


r 

i 
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<\  M  vrpv  uii  rvW/  prtATTA*^ 


Whi.e  it  adds  nothing  *o  the  present  forrul*t*on  of  •«  equations  cf 
ootion,  it  h  ^  ''nljable  check  on  the  basic  formulation  v>  investigate  t^e 
vorc  and  ec-r^y  relations.  Iking  (15),  (27),  (23),  (*»£»,  ard  (64)  ..o 


£-££**  ■**•#* 


-t  £  va^.  ^  •  (i/-n  *  $§*  v  \u$) 

-  V-F-t-il‘6  +•  N^WjCj* .  (99) 

ard^~  are  given  by  (34)  and  (85),  and,  If  W  1*  vorfc  done  by  the  external 
forces, 

=  J|ni(v%n«sf<*^*v^  o1) 

=  7-r+ri-(j+Qo+Qj^’;  deo) 


use  having  been  ^ade  of  (47), 

By  substituting  tbs  suffix  O  for  j  in  Equations  (4j)  thru  (49),  Jt  find 

that 


N, + 0,  +P.  =  Q„  .  (wi) 

Because  of  this  and  (49) 

N„+  0o+-P. ■vfNj+Oj+P})^  -  Qo^-Qj^,  (102) 


f 


V  F-a  G  •*-  W.  +  O0  vp<j4-(«J¥0,-c?J)6|'i 
=  VF+S-G  +  Qe+Q.q* 


ato  (103)  nee  (&*j,  ( 6 a),  ana  /cants  la 


dT  j,  4'J  .  <iV ..  d.W 

«(i.  ‘  it  ejt  3R;  > 


(lea) 


T  +  U+V  =  W-*  C0MS1AMT,  {10« 

vnich  sirply  states  the  Tact  that  the  orfc  don*  by  vh-  exttraal  forces  Bust 
be  either  stored  ir.  the  fora  of  ciuetlc  or  potential  eners/  or  dissipated. 


6.  PRACTICAL  EXPRESS  10!)  OF  THE  INERTIAL  FORMJL". 


Not  all  _f  the  foregoing  equations  are  needed  in  the  computational  phase 
cf  investigating  the  elastic  deformations  of  a  vehicle  ir.  flight,  but  those 
that  are  essential  for  this  purpose  (in  Sections  3  -and  h so  far  having  been 
only  rath***  «bstraetly  expressed,  need  to  be  presented  in  teres  that  are 
suiwao  -  for  practical  use.  Among  thcr.c  essential  equations  are  (53 ),  (?•'), 
(ch)  (go\  r--i  (£0),  and  they  are  full  of  partial  derivatives  of  (or 

its  oonpcneut-i)  with  respect  to  the  generalized  coordinates.  For  tne  purpose 
of  conrputation,  these  partial  derivatives  need  to  be  expressed  in  detail .  For 
convenience.  the  subscripts  t  and  -A.  are  temporarily  dropped,  and  the  basic 
notions  of  equations  (1;  thru  ( 1 0 )  arc  developed  and  extended. 

Let  us  introduce  the  vectors  h0  ,  ,  5?  ,  and  havini’  co  do  with  the 

linear  and  angular  velocities  of  the -%jycoord£nate  system  relai  -e  to  the 
system  and  defined  as  follows:  * 


b 

If,.- 


I 


ty-.. 


u 

f- 

■#4' 

ifv 


(106) 


(lov ) 


(108) 

■jti 

(109) 


Tip  and  Ti,,  are  the  partial  linear  velocities  with  respect  to  <^°  and  of  the 
origin  of  the  Jjf  coordinate  system  relative  to  the  X-  ayatea.  CST  and 
SRj.are  the  partial  angular  velocities  with  respect  to  cP  and  o*o?  the  ^4. 
coordinate  lystem  relative  to  the  system.  ’  ' 

There  should  be  no  particular  difficulty  in  regard  to  the  linear  velo¬ 
cities,  but  some  diucussion  of  the  angular  velocities  is  definitely  needed. 
First,  let  us  note  that 


Jr  *  Os  -  Srs  “  ^  W*'*n  r*S 

-  0  wHat\  Vi^s 


Then 


(no) 
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HS*?*'.  v  ¥~?"  v *&*-*?" 


Also,  lec  us  denote  the  components  in  the  4r  system  cf 
t  ’ '’>'  1  y  a:;  u-vl  £>/£.  Then,  from  (103)  and"  (110), 


C?<o  end  o<-^  respec- 


and  likewise,  from  (109), 


These  equations  enable  us  to  write 


0<o  =  'lC°<,o  AW4  =.  ol'x 

(113) 

1 

it! 

and  to  show  that 

^0*i;  =  ^  an4  %Tr  =  ^ 

(114) 

ifl 

In  further  anticipation  of  terms  to  arise  in  the  equations  of  motion,  we 
derive  from  (ill),  (112),  (113),  and  (ll4)  the  following  relations: 


—  CX  ^  o  OC  £  .  (U?) 

By  symmetry,  this  can  be  extended  and  generalised  to 

*  ’&*.  .  (U6) 


Prod  (113)  and  (llU),  It  is  readily  round  that 

-  ft 

+O<ytxoici  . 

F.li«,.ur..:on  of  ^.0*  05(k  between  (116)  and  (117)  results  In 

*0  -  J' 


Sq* 


Likewise, 


w>> 


and  substitution  from  (118)  and  (119)  *nto  results  in 


&-fse  =  ^ 


(117) 

(118) 

(119) 

(120) 


The  derivation  of  (115)  thru  (120)  was  of 
of  them,  the  suffix  aero  could  be  replaced  uy  a  latter.  In  the 
realm,  this  means  that  the  relations  expreseed  ty  these 

f-ahle  between  degrees  of  freedom  involving  only  wsail  motions  (type  (3))  ** 

S,  u  S“  U  .«»»  or  tn*  (2)  «*  ««-  «“«“  °r  ’■w  (3)- 

B.  «corO«K.  Vito  (6),  O  •  tot  iaif''  or  may  nor  b« 


dq*' 


aero.  let  us  introduce 


*V 


<rr 


(121) 


Then,  since  — ,y 

a5  »  ^ 


(122) 


«*  %-  *  Vr 

=  J \  *  iV-»'v 

=  4*1  *^r- 


(123) 
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*!??,.  -.0  '123)  facilitate  th*  writing  of  the  ‘‘oilcvi*^ 

A'fc  =  +42. 

Too 

_  ...  +■  0*  /  "v7  .  (12- ! 

iT  —  ii?  a. 

+  Ef 

—  Ky-t"  <J£  <125^ 

I-ct  iu  introduce  the  two  assumptions 

fj?  =  °  <xhd  X.  0  (US) 

firvl  exnre«*ions  <v»r  ♦*»  cccc-rf  p?rM*l  icri/atives  of  with  respect 
to  the  g-^erallzed  coordinates.  It  can  be  found  without  to?  ruch  trouble 
that,  under  these  assumption, 

^  -t- cta  x  ( 5^=0 )  («-) 

4^  =  »  -  |(128) 

*<fv  +  o*-v  *  0*.  +  cj<-**  t^i.  X’O)  V'J29) 

The  first  isportant  question  that  rise*  now  is,  hov  do  ve  equat-  {  ) 

with  (125)?  In  answer,  let  us  observe  tnat  (125)  is  the  general  fors  of 
expression  for  :  that  theZ^ni^  ere  given  as  possible  or  tenta¬ 
tive  rali-s  of  $pA>/&£iK  #  nna  that,  therefor  .  tfce$K^  ^Ul  be  given 

In  the  same  fore  aa  (12$).  Thus,  let  the  given  **r*car  velocity 
of  the  center  of  »cs*  of  section  (.  ,  and  let  its  given  angular  velo¬ 
city  be  £.  .  Then  * 

!.  <W  =  ?Kl  +  ff£U  +  5  (130) 

and  substitution  i ro=  this  into  (  40  )  results  in 

+3 lu. -V Jrk 

=  &tii  +  W  +■  +-’rr\.<.')  +■  <Tv.t<+^(L 

=  h*;  +■  .  «■  .v.fc  (131) 
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“c  ‘  nM  -  r  *  O'.  U  *  * 

+  '^K  _0-33) 

.«_ .  «s  —  /  consider  the  practical  problem  of  oc.l.,ng  (^3 )  ror  - _ ~r  . 

«  *  *he  use  of  for  4.  ana  of  Pq.  ..on  ( 130),  the  first 

terss  "  itqj  bee  oses  U 

^  tt  t  ftU ( §i  -t-'GU)  ^ ( Ki  -•-  i  *  tU) 

=  £.[fci$  x^icir  0(x£.y»i.*  dr^ 

p  p 

f  *  {&;. * t  m^TJeO  -  ^ 
t 1,  K*  ^ '  ( <**.*.-*-  K,  *  53  **.)]  .  (13* ) 

Pros  (2)f  (2U),  and  (121),  the  *olloviftg  results  arc  ooteoned 

=  O  ,  (135) 

P  ~v  P‘ 

=  **  u%6) 


This  :i.isinates  three  terra  of  (ijt),  axyj  5ta  final  tern  is  transferred  as 
follows 


=  J.*u*  [&. (TIvV 


.30. 


£<snc»n* 


-  j;,  b'Z  £m(U  Sr,  v,v  \>;\  *ui\-  Mf;.  •utV') 


T/  o 'S  / 

JM  ^  *C  , 


U.V7) 


wi,.,,  •  ’rti =  |.Yn,<(6r5'u'iX  ^  -  *w<5.  ■Uil') .  <*#> 

Tut'  are  easily  cesn  to  be  the  r  orients  anu  tue  negatives  of  -‘.lie 

products  of  inertia  of  nee  "do -.  i  about  its  out-  uxi  - 

For  convenience  in  treating  the  remaining  tcims  of  (l?,*.*,  we  introduce 
the  fenautaticn  symbol  *  0  if  two  suffixes  are  the  aa-ue 

■  1  if  (rst)  is  an  even  permutation  of  (123) 

«-  if  (rat)  is  an  odd  permutation  of  (123), 
the  even  permutations  of  (12*)  being  (i.23),  (2_>1),  and  (312),  and  the  odd 
permutations  being  ( 2 7l 3 ) ,  and  (132).  Thus 


tftu  -  ^hS-t  ^(VT\C  &i  jf'x.L  f 

|('m^  Xki.  x  =  c™  Yu  y\r£  f 


where  A'3  *  . 

Substitution  from  (135)  thru  (l^l)  into  (13I) 
results  in 

5  %,  ~  5  (c»-»tTtV(\i  ®-t  fa 

+  C'r5t  *1 

=  Jr  I  (CwtlRue!^  ”t*  c$tw  ea;\“ 
+  e«(  Wit  a;c  j . 


(135) 


(1U0) 


(1*H) 


(1^2) 


(1^3) 
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7  •  :»  i'.;  -  -f  (<•;'  is 

= x-  A*  1 1  ^  &*  A  lfc“  &  f  ^ 

=  J--^  Irs  ,  (1U) 

where 

iri  *  li  (<>>•*  {h-t-  (M.^.,0 
-  |£inu  [u«:^;  <0K + et.  ^'1) 

-  (o-r  +  eC.  ^’1<')(©-!  +  a**  -u'&)] 

=|  («r®^e;ceU 

t-e;;  el.  |’,m^(&^vlV  tjiI-  ^'*)1 

"I,  Knots'- «-c®-’) <-<teV.HV.-J  .  (11*5) 

The  I,,  ere  toe  rooents  and  tne  negatives  of  the  product*  of  inertia  of 
the  vehicle  about  it*  axes. 

The  third  tera  cf  (1*3)  contain* 

a.*.*  =  | ^  »u  (^u-t-sT.ir  B,.  *  ■o.-O 

-  fj  WW  |V 

-  Tr  t  ,  a«) 
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jf-*)*  •nlMj* -ft 


VS 
k  IT* 


-  ^  V  ^  k.  (,<>rs  ^  ^ 

S^etiwition  frow  (lt»3)  thru  (1U7)  inu,  (43)  result*  in 

^£®<^c  +  c»t«<iAXKi  -11  ^  H'itl  j3*c) 

’*"^r  "^Kv  lf|  _  'j'i  *  Jt  f.  ft!i  pf^ii 

(fc*  if  f 


{1^7  ) 


W  o 


Noting  that  u,. 

the  fora  "  ' 


(1-6) 

ve  can  pat  this  in 


-f  A'iJ  +  e^ni^  B^t)  .  (llt9) 

U<  .  K,  .t 


L« 

>ntl  i "•Ik*S&; 

Uten  (ity)  becceea 


lV»"  (6  rsfal-  --  Cw  ^  Ht 

-cntL‘i.- E|(c*a*eJtA'£t +  ef{  , 


(150) 

(151) 


(152 


”“'h  5*?,^  *olv'i  foi  the  J\  by  feaill.r  technique..  With  the  eld  of 
(H),  (lk6),  (150),  end  toe  fed  thot  X-k  *5c-  Te CM.  a.*  ,  It 
It  easily  teen  that  *  ~  St  ■ 


'  ^-rst  —  JJJ  H* 


t 

i 

i 


-3)- 


tf  >  mmt* 


f  -  , 


+  >ls<.  t-  it.  *  £i\  'VJ'’i. 

- irl  -/*, 3*1V>  •*  ^v. r- it*,  s C  '\>'~- t 
=  iA^  ->•  e^La-'X  v  :^  3',\  y>:\)\  , 

^  ^SC  (.  ^""i<  --Vv  C-iU  3  jti  -*  ^  _  (1>*) 


Similarly  (1*0)  is  transferred  to 


T,  =  Jv-arKl<  +  *s  5  Jt  2-t  jft*.  +  TrC^ 

=  M^Vvt  +  c,.*  z-t  £*.  +  c;) 

■*"  es£  (  <J'*(.i»+£st<i.  $*,•  "W-O.  (155) 

•fe  nov  turn  oji'  attention  to  *he  evaluation  of  the  P  •  »  f 

defined  ir  (53)- 


?,!>•<■  f<~  $r‘  »  |i*«  (K<- M*W) 


+^"1;  +  e‘tl  (  OJ'l  +  ctu,  "Oct)]  .  (156) 


i 


(15  ) 


f- m- *  yu  c,„  =  c„h  g  , 

G-  =  £,£,  Yru  'r-< 


’'^hiIvv/^  “  iru  j 

= 1 &  *  <*  ^  e*  <n& 

=  £  e:;et,  f  vn^xi'^  <r4£l 


-  I  e^eL  A'r*  •, 


^J<4\ 


=  Ct 


tuv  I  e‘:  £jr  I  Vric^  (®-c  -*-  e^,-  -v'fx) 

.-I  ^.t 


(158) 


(159) 


(1601 


(161) 


=  C^S,  zlfijtfy  ?w  , 


(162) 


3 
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'■WOW’**.. 


r^-/.  = )  \>'£ 

»  X-/ 


6^-v  "  K/*v< 


-tbtlt-t!  O'  free  (»>')  t^nj  (*.6c’)  i  r'*„  ( )  r**6  its  lr 


?  £ 


=  +  C~i  A)  ,0>u  d-J  ! 

+  t  f  e;,  et\  A'/  +  Ct.  ►  e‘;  e^ 


Tree  tils,  «>ttl*v  5  e q iv. i  to  r  {&r«d  sumLng  over  f  ), 


y  &  "J 
£'. 


(165) 


+  /jw^r^l  +  e^etiA^V  -  er.eV  <£J.  /"£*■,  j 

-  vn^.:  c[  + 1  ( wi.&r  ^j{  -v  AjO , 

,lncc  Gk=6»--,  e;;efr£«i  <xmJ 

/>V,-  =■  /“irptC  . 

/^sj  lc  nov  ruined  fc>  cubatltntlon  frax  (16b)  and  (165)  into  (51)* 


-%■ 


p-v ;  ' ■" ^  CI " £,m‘ t a“j - c- 

-r  / 1  r.  .f ,.«►  /v"*' 

•  -  -j  -■-  t-i v  *•“■  ’  c->c-*'  -  V 


■'*  CtwrSu^itX/jt  fg3.1r>  ) 


*  ^Uysflq  -  cj)  -  C«,^6m. 

+  E([fn.  «  ;  -  °Vli)  +  5vs  A'ji-  e*.  A^c 

“’  Ct«*v-  ej(.  r^-iri  1 

*  M  i rs  C-j  ”  (3^c  C  j )  ~  c  $4 a.  -6-j  & hu 

+s*Lf;  -  l?  +ti[e;i ei (StuA'jT-  A'jc) 

—  Ctu,v  C^i  etc  np%n  ]  (: 

ut  ,  C 


and  substitute  free  (155)  into  (60). 

Ttec,  after  *ir.f lificsti;..  aaUng  use  of  (22),  (23),  (25),  (lj5),  (13^)> 
(138),  (1^2),  (1*»5),  (150),  (151),  awl  (159)7  rreult  It  obtained  tha* 

k,»=  m[cjc;+c„,(ci^-rC^')^*] 

+  /<jK+-/-]  /■!  Ir;VCjH*-rC*  Hj 
+  C(-St  (AT  L'k  +•  l”  )  +  I  ^ 4’. 


*  jsqaSftiiie.'o 


♦ Cmft;  £  e;.  a':: + ■£.;  i,  e;k  a',s:  ) 

a“) 

1  -^-j  I(  £5,.  B  ,.t  VWtc  t-A-U.  ^  ft  '*  u 


+  A  ^ x;  Wksi 
W 


(168) 


TMs  equation  *.s  coxfirred  by  rfritln*'  ojt  fi'J  ex;re*slo-  for  ♦he  <ireti«*  »*nerr, 
x  and  uxlng  uae  of  \72). 

To  provide  icrc  compact  notation,  let 


PjK  = 

%,*n'lTs'-  If*'  > 

(169  > 

Ak  ’ 

5  Cs»*  £  el.  A*v  i 

ul 

(WO) 

= c^s-.c:  +/j;:  a”  ); 

(171) 

w; 

*»<« 

a7?: 

H|K 

h;,{ 

<173) 

TKch 

m,k  ---  mfcrc; ■*-<:,.«  ] 

+  A*  \  \-z  V"  C  j  Vi*  -vC*  YA^ 

•*■  (^Kt  L*i.  +■  A.fc  +■ 

+  ^(c,5a‘;-A-j  +MS) 

+  PjK+‘Aj|<.+  Wj*;  A-  -«.jK  .  (171-) 
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If  vw  nuoswriut  i..  (12C)  Is  rerln-ed  by  L  and  the  res.  lt  applied  to 
( 11-'?},  !  t  'an  be  demonstrated  that  interchanging  K  and  L  does  rot  change  the 
value  cf  (129',  vh'-'i:  is  ns  it  ohould  be.  iv.is  "symmetry"  cf  (129)  depends 
■>!',  the  retention  o?  tne  last  term;  but  X-u  .  thich  * .  n  factor  in 

the  ’art  term,  I  :  niffieult  to  obtain ‘and  likily  +o  be  3i?all;  tne rc fore ,  a 
means  of  droi-dn*  5 1  tut  without  destroying  the  aynnetry  of  t'-->  e-tuat1on  is 
sought.  T.i.v  Is  accomplished  by  a  simple  “"eeegir.g,  as  f town : 


0^  '  ’  t (J^L  +  iur'‘ 

=;  c*K  x  &Z  +-^v.  *  C*  +-  x  (otu  x  aj)  c\t  x  (c  - ,<  x  a})] 

x^- (o^ ^  + ('6^53)5^ 

Bnploying  (125)  and  (175)  in  (59)  results  ir.  - 

_  N  Pi  _  _  %  f  _ 

ia^l  -  £.  ^  Vni<(kji  +  0JU.  x  -utco  *  !<**<.  *  (176) 

. .  'ij 

If  this  is  now  expanded,  the  result  is 


where 


Ktdfe  t  (tv*  (C  3;"  +  «:v  s;‘:)+ <t  ocLt  d';t* 

■*■ (b^V  <f  ii  />-'v5)  o*fiO  , 

S/*?  = cjiiT  , 

--  £  ,irrut  ( e,s  v&  qci  -  -xl  <rrr< ) 


(177) 

(178) 


-  f  A'**  -  A" VS 

""  <->  hs  jo  i  X  jv  ) 


(179) 
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✓  rS 

«!Ki  U- 


"I  L0*1  J  L  (  *"  St*  ^  »"»•*  <rfc«A  ) 


-r^-rTi 


(iso) 


xi  ve  note  that  ^jK;-  *  0  when  j  nri/c,r  K  equr-i  zero  t»no  that  jf3  jt  s  O 
when  i  ■=  O  >  then  it  follows  frcm  (1'C?)  f'i,a  f- 

U 


OKj 


=  ?  (c. 


X-  v'-hSi'^0t 

i«l 


•VoC^,.  oCo£ 


a'") 


(181) 


and  6gTN-  J  oCiO^o*  <X'fl  . 


(102) 
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.  «t4W«  *»' 


PRACTICAL  EXPRESSION  OF  TILE  AERODYNAMIC  FORMULAS 


The  aerodynamic  forces,  being  external,  are  accounted  for  by  the  use  of 
Equation  (47).  In  this  use  of  (47),  however,  only  tnose  particles  lying  on 
the  curfcc.-  the  vehicle  will  be  involved.  In  the  present  formulation,  we 
have  re«-t»ui;-e  to  the  simplest  available  »«rody..aaic  theory  tiiat  offera  suffi¬ 
cient  rcr- •(  i  -j,  namely,  Newtonian  flow  theory.  Let  R  be  it  unit  vector  located 
at  a  certain  point  on  the  surface,  perpendicular  to  the  surface  at  that  point, 
p.nd  pointing  outward.  Th;  velocity  of  that  point  it  &  ,  cut  given  in  (15; . 
Lotting  b»  the  atmospheric  density,  than,  according  to  Nvft  'nl.vn  flow  theory, 
the  aerodynamic  force  per-unit  area  (7)  at  the  given  rvint  >n  tie  surface  of 
tut  vehicle  is  as  fellows: 


1.  When  r\  •  if  — 

2.  When  H  >  if  > 


F  =  O. 


r-  -  *e(R  v-)3- . 


(183) 

(184) 


The  scalar  Tv  V"  may  be  called  the  "piston  speed"  of  the  given  point  (or 
the  down  wash  at  that  point)  and  Is  symbolized  by-tf/f  thus,  when  ■ua'^q, 

T  »  -It  £  2..  (185) 

For  the  purpose  cf  evaluating^.  It  Is  noted  that  the  first  and  last  terms  of 
&  as  given  in  (15)  are,  under  normal  conditions,  such  more  significant  thun 
the  two  middle  terms.  Dropping  these  less  significant  terwi  results  In 


ar"R 


(186) 


A  considerable  practical  advantage  can  be  realized  if  -Op* in  (185)  is  replaced 
by  a  linear  approximation  (or  expansion)  about  the  elastically  undeformed 
configuration.  Regarding  fir  as  a  function  of  end  6*  for  this  purpose,  and 
using  the  subscript  O  to  denote  the  undeformed  configuration,  we  obtain 


«r*=  v*  +■  2.  -ur.(^  *|«  +  ^  Y)  . 
TV.  »  V*  K  i 


and,  therefore. 


(187) 


(188) 


tu-*=  (v-hf  +  av-'R(7-^  c^4-n  .4*.  c^y  (189) 
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In  applying  [‘•l )  to  the  calculation  of  the  genei-niized  as*  1/nas^c  for-c 
ve  are  ispe*ieu  to  replace  tn<»  surae^ion  ever  ^  with  or.  \r*  2,{retion  over  tn* 
surface  of  t*»e  l-th  section  by  virtue  of  the  fact  that  on.y  pilots  on  the 
surface  a*e  involved.  Lett,  be  the  surface  the  i-tn  c^-tlon  then 


=  -  l  (  ft  •  h  ?  -ur  s 

_  r,  f  *  r - ,z 

"e 


+  (W) 


where 


d?i) 


The  following  developoent  of  formulas  serves  to  eaXr  this  sore  practical 
for  numerical  cooptations 


h  =  J  ■;  n-". 


(192) 


,&  n’r+Vr& 


4n  _  cJj- 

3f  E|* 


AS' 

:  *  x  h  -V  It-  KT 


-  V'  ^ r 

-  iy.  l/.*  f\  ■)•  _ 


(193) 


;n:,rtl„s  tnia  Into  i.,-A  *«  tr-.n forlne  1.  In  other  »«,»  »««*•  ** 


(J.  n  -  e  r  f  £. Tv ‘  Vs (in- Jit's (4t  ’  '*"')  n' ' 

*■  v»;  »  ^  ** 

i  2.  'I'V^ljy-TU)  (  js  T'ncV  ‘  (Ct'-H  c''*  V\-  *+-  )of 


*.  <»* .  f  >•*  wt’^a  ^ 


uYier* 


(195) 


-  tv;v:  i  e;, e;;  ]  j, 

i.*l  'Sl 

-2.?vv£  e;,  etc  (c*v* '*'*’<  (  n^V**!**® 

V  -  I 

-a?v”I  e*<.  f  .  n*  &  i„  <<*  qK 

i*|  ' 

-  -((v'v1Aj‘+J'l''',‘®i*cl'l+iV''<:'i|lH11)  '  (1?4) 

er«-t«iet(cw  <<£ 

+ln/jv'/t  <**)  > 

-  c;  =  £  ei  kfc  As , 

Substitution  from  (125)  into  (19D  results  in 

|r  vi  •  (itj *  ^ 

=  Y\'*(eSWJ  4-o~j  *'+  c  •f-3'  ^ *)  • 


(196) 

(197) 


(198) 
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I'.  !<•.  (or  .« .-.trof.ucc**.  .  a  .  1  i  >  on  tlv.  ■  -  .Mcc/'ilit  of 

tM  ’  '  i:  •;  tv.it  all  V'  .-Ot  o ,*cc  .uv  fi recto *.  V-.  thro-t  "vccte.-i  v  ” 

(that  '1.  ..cting)  notsles  that  _'i  .mch  nozzle-’  :;y>.**Vrlo  with 
rci'jw.1.  '  Is"  axis  (or  lire)  of  thru. L .  'et  erding  .ueh  t’.irur.t 

'.c  .,iu  wit  movable  structural  recti on  luent.lfiauj.3  \iy  a  subscript 

1,  we  r>3 ur.c  therein  a  tr?m1  of  wit  vectors  J*„,  with  ,Tj  poi  \ti  .r  in  the 

direction  of  the  thrust  and  coincidini,  with  the  line  of  thnnsU.  Jiscaun  j 
cf  the  sjvjxn.tr:'  of  the  nozzle,  its  center  of  ie.os  wit'  lie  j..  ub.e  thrust  \xis 
and  the  usual  conditicn  tliau  the  urit.ii1  of  the  triad  he  at  the  carter  of 
oass  of  the  section  can  l*  car  olled  •  ith. 

J&L  =  =  K . 

At  tlie  center  of  mass  of  a  section,  d^l  •3<y  * 

Representation  of  the  thrust  force  at  the  1-th  nozzle  as  J’ .  V  and  substitu¬ 
tion  Into  (^7)  results  in  the  fo'lovlnc  expres  lion  for  the  (generalized  forcer, 
associated  with  the  thrust  forces: 


(§t  :  1,  Kji  *  ’ili  fl 

=  £  'iv-i*  J|L  VvjcTi. 


(199) 


--7.,eTchVTc, 

E  beinc  the  nur.be r  of  thrust  vectoring  noztleo  (or  "engines"). 

k  A3  has  already  been  indicated  in  (7),  the  ei.  are  functions  of  ^  and  the 
9  •  The  correct  Inclusion  of  the  dependence  of  the  e!h  on  the  9*  would  be  the 
best  procedure  and  would  enable  the  program  to  reveal  the  interaction  between 
thrust  and  elastic  deformation  even  to  ths  point  of  detecting  instabilities  l-** 
anj'  existed  .  However,  doing  this  would  iiapose  considerable  additional 
difficulty  and  go  beyond  the  scope  of  the  program;  therefore,  ths  dependence 
of  the  e*^  on  the  <)k  will  be  disrega.?ded  be*c.  (to  toe  other  hand,  their 

dependence  on  Cf°  must  be  and  is  included  as  shown  in  the  following  section. 
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gr^CTior;  cocrrT::.;  cr  tcvakle  st  sUC^riAi.  skciiq;; 

In  addition  to  thrust  vectoring  nozzles,  there  are  such  movable  structural 
sections.  a6  control  surfaces  of  various  types.  For  the  sake  of  simplicity  it 
is  assumed  that  the  large  motions  of  all  control  surfaces-  consist  of  nothing 
acre  than  a  rotation  about  a  fixed  axis.  It  is  convenient  t.o  place  tne  7^ 
veet*>i  oj  such  e.  section  parallel  to,  but  not  necessu  j.ly  on,  this  axis 
of  rotation.  Doing  this  makes  it  possible  to  employ  Euler*  *n  angles  to  define 
ttu  or .  -  *\^ion  of  both  thrust  vectoring  nozzles  and  control  surfaces. 


These  -.  ..gles  are  shown  in  figure  1  and  defined  ac  follcva. 

a).  *■  angl  of  rot  at  ice  oi  the  plane  and  axis  of  sw:.vc:‘  about  the 
r‘  y»  axis  (JL). 

A^1"  angle  -f  svivei  of  nozzle  (or  the  7*  vect..r)  ahmt  an  axis 
(l3)  perpendicular  to' 7.  and  making  an  angle  d)^  with  7-. 

5i“  angle  of  rotation  of  7'  and  d*  about  7* 

cl  1 

Ifcr  familiar  processes  of  vector  analysis  and  classical  mechanics,  it  is  known 
that  the  direction  cosines  relating  the  7*  vectors  to  the  7a  vectors  are  the 
following  e®^:  r 


<- 

cosX( 

e*r 

-sin  Ai  cos 

«ln\j8lnS^ 

e*- 

cos<X  sin  Ac 

coa${  cosXi  c°86i  " 

sin  sin  Si, 

-cos^j  cosXv  sln&  - 

sin  ^  cos  Si 

e,r 

sin^i,  sin  At 

eir 

sin^coa  Xjco *Si  * 

cos^i  sinfc>£ 

eic 

-cinfi  cosXi  sinSi 

+  cos coo  Si. 

(200) 


It  is  understood  that  Ai,  and  Si03™  functions  of  cf  as  determined  by  auto¬ 
pilot  or  flight  progransaer  connands.  ' 


i 


DO»K 

HOfCt 

nj'lj  AHOT,  AM  l  .  ROT.'.Tt  ABOUT  AMO  A8C  l  TO  7. 
k)  J's  AJO  A  A?  i  fc-  -.^ayiT  A»Pi«  *  TOTj 
«|  J*,  AMO  ij,  A*t  1  .  ftOYATt  APOUT  AMO  AM  *  TO  7,', 


figur#  3.  Orientation  AngUi  of  Mevabl®  Section* 


FORMULAS  FOR  THE  STRUCTURAL  LOADS 


Thd  shear  force  S  at.  a  specified  location  or.  the  vehicle  ia  ihc  nogo- 
ive  of  the  fium  of  the  Internal  forces  exerted  by  all  tne  part  Icier  of 
the  vehicle  on  the  particles  located  on  one  side  of  i  he  chosen  shear  pl*r.o . 
For  the  ttt»n.e  of  economy,  the  side  of  the  shaar  plane  -  ieeted  for  this 
purport  Jill  l—  the  aide  on  which  the  smaller  number  of  particle*  ie  fo.n»d. 
7nii  «  .  ntuAlly  be  the  side  away  from  (or  lying  outboard  of)  the  center 
of  mass.  Lev  the  absence  of  specific  designate  ->ti  as  to  which  particles 
t*.ai  section*  nre  included  in  a  summation  be  understood  to  mean  summation 
over  the  phrticlet  on  the  chosen  aide  of  the  shear  plane.  fLe n ,  with  the 
aid  of  (44) , 


I£LLF-*4.  =  F.X:  (&-m* 

c  A  *,/  *'l  °  1  < 


(201) 


Except  for  the  number  of  particles  included  in  the  susnation,  the  last 
term  of  (201)  is  the  same  as  the  right  side  of  (27) •  A  practical  sy.'snet- 
ric  expression  for  is  given  in  (175)  •  This  can  be  used  in 

(27),  which  in  turn  Is  to  be  used  in  (201). 


The  bending  moment  ^  at  the  specified  location  is  the  negative  of 
the  sum  of  the  moments  about  a  point  "2.  in  the  shear  plane  due  to  the 
Internal  forces  exerted  by  all  the  particles  of  the  vehicle  on  the 
particles  located  on  one  side  _of  the  shear  plane.  In  like  ox  inner  to 
that  employed  in  determining  £  ,  it  is  found  that 

V  *  1  *»lpl  * 

--ix- 

C  "F- 

=  L  J  Fl- l  ~  £  jf l*  x 

-pS  .  (202) 
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,:ir 

Li  i 


.  i 


CopV 


Fxrnj."  Tor  ‘ ' ,  -.tent  o!’  th*‘  -unnn*  ion,  the  next  to  the  lua*  ter  r.  of 

.  is>  ;  j tu*;  l  i.^iit  cicic  of  (*  ',  which  car.  he  used  ,-j  -xpnivi  this 
tert:  i'o»’  ;:i'Xcttenl  use. 

^>  -  riy-Sr  Ohd  TA  =  >i  y  V\  y  .  (203) 


The  ee.u^i.  numerical  uunntltiea  to  be  computed  (ire  the  components  $y» 
of  J  ana  Mj.of  It  is  uoSuu«.'o  for  th.  -ur poses  of  this  nronrj*" 

thnt  utie  selected  shear  plane  will  be  perjendlcular  to  one  of  the  J. 
vectors.  The  choice  of  the  sheer  plane  will  effect  th-:  interpretation 
of  the  results  S».  and  (  V-=  IjZ  •  Tr'~  ••w-epte.  If  the  shear 

plane  is  perpenalcuiar  to  T  ,  then  1  la  the  •.o»i.-o*.^t  .  of  the  shear 
force  In  the  direction  of  th?  u}  axis,5  g(  le  the  ahai**  force  In  the 
direction  of  the  axis,  5  x  la  the  normal  force  (being  oerpendicular 

to  the  shear  pi  anal,  if  the  bending  do  sent  about  an  axis  parallel 

to  V  >■  »A,  la  t,h<»  bending  moment  about  an  axle  parallel  to  a*  , 
and  the  torque.  w 

It  Is  a  prerequisite  to  (27)  that  the  sum  of  the  internal  forces  exerted 
by  and  on  all  the  particles  of  the  vehicle  equals  zero.  It  la  likewise 
prerequisite  to  (26)  that  the  sum  of  the  moments  about  any  point  In 
the  vehicle  due  to  the  Internal  forces  exerted  tor  and  on  all  tbs  particles 
of  the  vehicle  equals  zero.  These  facts  are  deduced  from  Nsvtoo's  third 
lav  of  motion;  and  It  follows  from  these  and  the  definitions  of  S  and  M, 
leading  to  (201)  and  (202),  that  determining  3  and  M  bf  sunaring  over  the 
opposite  side  of  tbs  shear  plane  should  change  their  signs  tut  not  their 
magnitudes. 

It  has  been  noted  in  Section  3  that  (27)  and  (26)  are  not  satisfied  In 
the  SLP  because  of  the  assumption  that  the  elastic  deformations  and  fuel 
sloshing  motions  havs_a  negligible  effect  on  the  large  motions  of  the 
vehicle  and  on  F  and  0.  Furthermore,  the  aerodynamic  theory  etqplcyed 
bare  (Newtonian  flow)  la  different  from  that  employed  In  the  S1V  program. 
This  difference  between  the  two  programs  further  Jeopardizes  the  agree¬ 
ment  between  them  as  to  ?  and  %  and,  hence,  the  satisfaction  of  (27) 
and  (26)  In  the  ST.P;  therefore,  it  canr”*  be  expected  that  summing  over 
the  opposite  side  of  the  sheer  plans  will  satisfy  tba  theoretical  require¬ 
ment  at  changing  only  the  signs  of  3  and  M.  This  represents  a  failure  to 
satisfy  Newton*  s  third  lav  of  motion  and  may  prove  to  be  a  rerlous  defect 
In  the  program. 


•  iromsioti  or  rjzL  sloshes  prscTb 


The  basic  formulations  for  fuel  sloshing  as  o'jtMloed  f~oa  available 
literature  ere  presented  in  Appendix  Z.  Here  the  concern  Is  *i>v  to  incor¬ 
porate  :hs  effe-ts  of  fuel  sashing  in  the  Structural  lords  rrograe. 

Vlv»w  each  ttni,  there  are  two  directions  of  siosning,  designatee 
9oe*vsp  _  '«eiy  as  longitudinal  and  lateral.  Vs  choose  toe  letter  u 
as  an  Indicator  of  tho  ulo thing  direction,  lor.,.tudical  sloshing  being 
’natcated  by  ‘tcting  u>i  ,  and  lateral  sloshing  being  indicated  by 
letting  u*5  • 


Tor  each  sloshing  direction,  there  are  tvo  pci.si  uodes. 

The  latter  S  is  the  node  indicator,  tbs  first  node  being  inflated  by 
lotting  s  *  I  ,  and  th%  second  code  csisg  Indicated  Vy  letting  a  *  . 

Since  there  ar*  .vo  sloshing  directions  and  two  possible  aodv*  for 
each  41  Faction,  there  are  four  possible  eloaning  degrees  of  freedom  for 
each  tan’s.  The  rusher  of  the  tank  is  Designated  by  the  iotter  » ,  and  the 
degree  of  freedom  la  designated  by  K.  The  following  fcrsmla  is  used  to 
deterslne  H  in  terse  of  i  ,  u,  and  S * 

k  "  4  (  i  •  i)  ♦  u.  ♦  a  (.«•->)  . 


The  following  tabulation  illustrate a  these  relations. 
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MM 

Mpt* 
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MM 

mam 

Mm 
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m 
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a 
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m 
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B 

B 

a 

BB 

B 

B 

s 

m 

E 

m 

*ft>a  prograo  allow c  for  a  oaxirua  of  ten  tanks;  therefore,  the  largest 
possible  valut  of  H  for  e  fuel  sloshing  »o**  la  *0.  The  nuatar  of  the 
clastic  degrees  of  froedoa  for  structural  daforaetioc,  control  surface 
*o  tat  loo,  and  so  forth,  starts  with  *1  and  proceeds  to  a  eaxiww  of  57, 
t  giving  a  possibility  of  17  "structural"  Cegraes  of  freedoc. 

The  greatest  problem  thot  arises  in  connection  with  tho  effects  of 
fuel  rloshing  In  the  Structural  Loads  Program  Is  the  corputation  of  the 
terae  MT*7Ti33),  AV(lW)f  end  SjJuUTS). 

In  Appendix  I,  fcrwulas  are  given  for  the  effective  ooeents  of 
Inertia  of  the  fu*l  tank  axes  for  rectangular  and  cylindrical 

tanka.  The  equivalence  between  these  and  the  H'T%i  is  as  follows,  the 
sUbecrlpt  P  denoting  fuel; 


*9* 


H 


'i 

f  '  3  .i,  '  I  ty  i 


o  wVien  r*s  ! 


for  a  rectai^^l  .  wans. 


H 


H  r»s* 
H  fro 


H  Vai.  « 

H'ffn 

wa  h;,„ 


■  J,«.  -is,. 

r  O  v/V\iv\  r  ^ 


Jri>  * 

-  0  V*4  -•  O 

*  o  when  r  *  & 

•  O  for  a  spherical  tank. 


*o*  a  horizontal 
cylindrlc’j.  tank. 


xor  &  vertical 
cylindrical  tank. 


(2C6) 


<207* 

(206) 


In  rectangular  tank,  and  for  longitydin.il  noshing  in  horizontal 
cylindrical  .anas,  a  spring-aass  oechanteal  analogy  it  uacd.  Each  oats 
in  this  analogy  haa  notion  in  one,  and  only  cat,  degree  of  freedos; 
therefore,  it  can  be  identified  by  the  subscript  k  ,  in  accoi  dance’ vitb 
equation  (20*).  I ike  vita,  the  location  of  ,r>*  it  given  by  the  coordinates 
*i*y»  »  *  *  •  Is  the  case  of  longitudinal  oscillations,  X  .  ■  q  *  and 

y„*o;  for  lateral  oscillations,  x**o  and  ;  in  either  case, 

£,is  sirpljg*  ,  a  constant. 

For  the  purpoeo  of  evaluating  the  A.<  and  the  5,*,  }  lt  is  necessary 
to  relate  the  raises  and  coordinates  Just  discussed  with  those  appealing 
in  <'il*2)  and  (173)*  An  inspection  of  (2T&)  quickly  discloses  that  the 
particular  east  particle  vithin  tank  »  is  idsntlflcd  by 


ii  *  u  ♦  al4*0,  (209) 

so  that 


K  •  *,  l  l-O  *  it  (210) 

With  this  relation  betveen  the  subscripts  »  ,  ,  and  k  established,  it 

is  clear  that 


H'l  ,X  ’  "1*  (2 11) 

‘  ^  K  wvVvcvi  u.  *  t  ^  (212) 

*  o  wktvv  u*aj 
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^  u 

y~  °  wV^-'1  a  i ) 

q'  «Ufn  u  3  ) 

(213; 

*•  t 

£  «  *  c  cias1c*v\t 

(215) 

Making  proper  applications 

nov  results  In  the  : Ursulas 

OV*  * 

|  whi^  u*i^ 

*  o  w»e**  .* 

(3iS! 

•* 

<^a.X  * 

d  y4  *  *  C  wV,«*v'  a*>  ' 

'  s  i  wUew  a  a  J 

<217) 

<r;v 

d  l  ./  dcf  ‘  *  O 

(213/ 

dusserlzlng  (212/  thru  (2161  ie»ulta  in 

v.;  • 

q*  when  r-u]  r  t  5 

0  vihe*i  r»u  1 

(219) 

■ 

H  „  v/VreiA  t  •  i 

(239) 

<C  • 

l  wV\en  t_  -  vx  ■) 

-  0  otiwwise.  J 

(221) 

Cubit i tut Ion  froa  (211) ,  {219),  (220),  and  (221)  late  (1^2) 
lotting  the  un'enovr.  q*  •  o  for  thla  purpose  results  la 


A  *  vri  *  2  .  wView  x  *  3  <xwJ  t  *  va. 

■  o  o4ViCrw»te.  (222) 

Slcil&r  substitution  Into  (178)  results  la 


•  m*  wWev\  r'S-u. 

c*«d  j*  v«. 

•  o  o-V^erwtfcC.  (223) 

*  Tbs  reader  Is  reminded  that  (222)  sad  (223)  are  applicable  only  to 

rectangular  tanks  sod  to  longitudinal  sloshing  In  horizontal  cylindrical 
tanks.  For  other  tanks,  the  A* L*  accused  to  be  non-existent,  and 
the  V4'*  ars  sore  or  leas  circuerented  by  arriving  ai  the  41)  K  (which 
equal  £  )  by  another  process. 
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s 


The  formulation  p’  escoted  xn  Appendix  I  for  lateral  oe^lHatlons  in 
horitxjntf1  cylindrical  tanxa,  for  vertical  cylindrical  tauks,  and  *"or 
spherical  tanks  lead  fco  expressions  for  the  kinetic  e-id  potential  energies 
of  fuel  slosuing  rati  nr  than  a  spring -*3a*s  Deccanic'*1  analogy.  Once  these 
for  the  kii-rtie  e^rgy  are  extended  to  e* count  for  the  other 
z*  *  >on*  of  the  vehicle  at  veil  is  « «*  sloshing  o*  the  fuel,  they  can  t>e 
used  a-  «•»  (72)  to  obtain  expressions  for  the  contribution  of  the  ft«*l 
to  t  •«  M,k. 


To  lateral  a loaning  In  a  horizontal  cylindrical  tank  ve  affix  ,he 
ruoscrlpt  •  toT,  M,,  ^  ,  and  a  und  saxe  the  fol loving  tebstitutionj 
in  the  final  equation  giv-* >  in  Appendix  I  for  tne  kinetic  eosr© 


%  5  V 

i  •  is.. 

A  >•  v  •  A  %, 

&)4l  *  Svt 

[""I' 


(221*) 


In  addition  to  this,  vc  introduce  a  transformation  cf  coordinates.  Let 


■fL  •  +  l  .  -  O  ‘  U. 
then  \,20k)  can  be  expressed 


i<  *  .jp.  >  a  (.s-i) 
Q 


(225) 

(226) 


S"  -  Jl* 

»  A.. 


..  (V-  ) 


Tnio  results  in 


T,  •  X, 


♦  tv,'*  X  Snqr 

Differentiation  of  T»  results  in 


r*  U-\) 


£Jl  .  ,3  a  n.-k)  R*. 

dqr.«t*-A 


v  c?  cli  i/v»  R**,  'll/  &«, 


(227) 

t  U-U 

(?29) 

1229) 


_ _ . 

d<\  t  '*tv  ° 


Oq  ^  ‘ 


a  P.  ”-.  5-.,  R*..  A  t.  /ai, 


.ae.a.i-.RU. 


j  oeooc2n^  o  "structural*  (that  la,  non-fuel  slrsm  ae^ree  of  freecLc, 
and  .,  being  de  xustrabiy  equal  to  — — ■  by  (62)  'i?4-),  - 

recognition  that  here  <j*  ‘in-i  ■\j  eq«a»  f~c\  9 

Frees  (230),  for  the  wse  lr.  which  J  anu  u  >  -  -  i  sloshing 

In  a  horizontal  cylindrical  tank,  ve  define 

u)k  •  a  e,  a.  ^1.  A  t,  /x  t, 

w^v\  j  •  k  •  C  s  -  \ 

-  o  wWevq  j  <  k  (; 

Froa  (23i),  for  the  case  In  which  W  denotes  lateral  sloshing  in  a  c:ri- 
zontal  cylindrical  tank,  ve  define 


’  S  f,  a*  lr„,  Rt.  Bt) 


For  sloshing  In  a  spherical  tank,  we  affix  the  subscript  *  to 
T  1  M,,  f  ,  a  ,  and  R  and  cake  the  following  substitutions  in  tlv.  final 
equation  given  in  Appendix  I  for  the  kinetic  energy: 

<*'  *  v,'r  (  r  •  >  «r  a  ) 

C  *•>  •  C  » 

D,.  ,  -  D,,  (2: 

[virrr  ]*  •  v*. 

In  addition  to  this,  ve  use  (225)  and  (226)  again  and  introduce  the 
follovlng  transformation  of  coordinator: 


^  •  v». 


ihls  results  in 


T,  -  i  M„(vr)**iire.a*  R^E  "}* 


R*v;rz'  D„  q 


1 
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-r* 

Dil  f.-rentiatiuu  cu  >,  results  in 


ST, 


pi  *  ^  I 

U  O  * 


-  11  <?'a1  R \-Tt-  >H  f 

Tre;  ^rOv, 


a  i.  i-  i 


(-37) 


_ — -  R'  c:*i  /  V  . 

SaT; 


a*  ^qt77^  *  nC-.  ^  R*  Dili,";  vi  yx  , 

j  denoting  e.  structure!  degree  of  i’reecom. 

From  ( -~,3 ) ,  fw*  tic  case  in  which  j  and  k  denote  elouhlna,  in  a 
spherical  tank,  we  define 


,oG\ 


(239) 


Wk  -  we-,  a1  R*  cst  / xv, 

whev\  j  *  k  *=  »■  a  Cs-\) 

1  o.  w  View,  j#k  (2Uo) 

From  (239),  for  the  case  in  which  h  denotes  sloshing  in  a  spherical 
tank,  we  define 


m'K  r!ds;  (sm) 

Making  uoe  of  the  nT*  from  either  (233)  or  (2Ul) ,  we  compute  for 
a  spherical,  tack  or  for  lateral  sloshing  in  a  horizontal  cylindrical  tana 


The  equivalence  of  this  to  (231)  or  (239)  should  he  noted. 


(2U2) 


Per  slushing  in  a  vertical  cylindrical  tank,  we  note  from  the 
given  equation  for  T  in  Appendix  I  that 

d'T  c 

<3^j>  vru*-  (243) 

L  j» 

Wow  y  and  %  need  to  he  related  to  the  coordinates  for  measuroment  of 
the  structural  deflections  and  the  sloshing  of  the  fluid  in  order  to 
determine  expressions  for  the  MjK.  For  this  purpose,  we  employ  Equation 
(2-16)  from  Reference  (7).  Putting  this  equation  into  the  terms  that 
are  appropriate  to  the  present  purpose  results  in 


•  !r'-r  fih  UaM 

Sr.  detersinlnf  the  partial  derivatives  cf  (21*4).  *«-  distinguish 
oetve-n  structural"  end  "fuel  slosh"  degrees  of  free1  as  before. 

Thus,  ;u.  structural  degrees  of  freedom  lassunlrw  u-£  ,  t*i*t  is,  that 
tv.«  >io  -  l«»  in  the  plane). 


a<f 

2)  q  i 


4*  -  Ji«, 

eut*  .  « 

dq  J 

ijr*-1  . 


For  foci  elo&cing  degrees  of  frecdoc 


cx  t  , 
ve^U{*{!  J 


dt. 

dq  r 

at1-,.,  -o , 

*  r 

when 

The  follovlng  substitutions  are  also  sade‘ 

til  ||  a  M  p,  t  ‘Vi  ,j  *  U  |K|  ,  W  t,4,3  .  A  tfi  (.S-l) 

^  *7  >  „  *  «  * 

•sa*1  0  ukt  >  •?!***  -u  ,  vnA  *  -u'i  ,  *  2  Cs-x'i  *, 


VTi  is  •  XL  H «  |  ,  Y>\  a 


e 


(8*5) 


(Sit) 


W %  u’-^.ats-U  ,  -e  is-O  (2i7) 

Finally,  substitution  froo  (3*5)  ♦  («46),  rod  (2^7 i  into  (244)  results  la 

*  Mfi  ^e,  rtj,  K„  h  ♦  u",»,  i 

*  cx'.',  .  u'it.C*  V.  *  « i«.VlK, 

4  -u'i,,  «*•  '«■  »  -U  J*.  )  (2W) 


ti  Cv-i)  ia,  Vlji  +  Alj  ,  JF  «i  l*"1*  otii 


*  aa  ,  ,  y.<a  Aa,  nj-,  +  aij.  ,  « 

*  u-'j  ,  -  4>„,  ,  (?<J) 
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•,**.>  1  --  il  Jt.ili,-.'  pjli.l  -  -U,,. 

■  O  wVlivl  j  *  • 

•  •'i*  mi  *iOsbX.  la  spfceri  »l  »nd  /ertlcftl  cylia4”l«^iX  tanfcs,  aad 
Ifl acre...  lag  la  horizontal  t&nfcs,  liiere  it  no  d;  oarlc  o&lar.cins. 


(2*0) 


’  2 .  POINT  0?  ROTATION  ASP  DYrtAXIC  BALASGIKS  CF 
CERTAIN  V?/  rpp  FLKaIBLZ  SSCT105S 


•ivM»  -very  section  is  *  ooint  a  (with  secilcacj.  po'it ton  v*c*  »* 
f  oi**  vehicle  jot  it  ion  vector  ?  w  which  say  or  nay  „c  coincide  with  the 
ccut.i  of  zazs  of  the  section.  Since-  dj  is  the  veoicio  rc^ltloo  ve«_U. 

Oi  f.  *»wr  of  cAst  of  section  •.  , 


X.  R.  ►  jl , 


(251! 


The  point  e, is  fixed  in  the  physical  saterial  oi  tv;  ».cv' »-  and  soves 
with  it  vhes  and  if  it  sc*e.«  Thus,  it  rety  correspond  >«  particle  of 
the  section. 


If  tse  »ecttc<«  ie  "oo^Galt",  that  Is,  has  cation  of  type  (£),  then 
7*^.8  i  point  of  rotation  or  the  section  -  a  peiat  in  the  saetica  that 
does  not  uoyo  relative  to  the  vehicle  but  about  which  the  section  ret  ate  5 
in  *  type  (2)  action.  If  the  section  rotates  ulout  <  fixed  axis,  ft, 
lie*  sccevfcere  on  this  axis. 


A*  for  actions  of  typ*  (3^,  the  unbalanced  notion  ;f  f*  relative  to 
tt*  vehicle  coordinates  Is  given  prior  to  that  of  any  other  point  in  the 
section,  it  the  degree  c t  free doe  defureu  the  section,  and  if  tho  section 
Is  eovabLe,  -f^idoes  not  save  relativ*  to  the  vehicle  la  that  degree  of 
freefioa.  On  the  other  hand,  if  the  section  la  "fixed",  or  if  the  degree 
of  free doe  docs  cot  defora  the  section,  t*y  govs  relative  to  the  vehicle 
coordinates  In  that  degree  of  frecdoe.  ’ 

7ro»  (150),  (135),  <13«),  Ml  (22), 


Pi 

T 

r, 


”*,1  o- 1*1  •  j~, ^ ‘ ^ "f  «*  " 


(252) 


t 


«■  -  -HT,  £. 


CT «l  S«  ~  -f  at  -  Pal  x 
Here  the  c^uikend  the  £ are  given  arbitrarily,  and  the  j  „  and  <T»ik 
are  determined  froa  thea.  This  la  necessary  for  the  satisfaction 
of  (136}  when  the  degree  of  frasdoa  (  k)  de fonts  ths  section  C 1) . 

-  _  _  _  — 

When  "Oix*  t*  *  *  than  t*jh*  x ,  cl  «i  *  ^«ik  *  J».  .  v  £  , 

and  we  Introduce  6 


X  *,  •  3  tu  * 


t" ' 


-f , 


(255) 
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When  the  degree  of  freedos  l*)  doe*  Dot  <>forx  the  section  O'  ,  * 

’  f  k*  •  *  (136)  is  still  satisfied,  rnd  x.  *,  »  ft.** 

Vhea  tse  section  la  eoYacle  and  la  deforced  oy  the  degree  cf  freedos, 
devotfilsc  the  a'I,K  and  tbc  {The  are  wtit’—ry. )  iV*s  cos* 

p*ite  ■'••■a  auoolt 


t 


Kl 


»n.K  <x  «»u 


•  °-'«,V-  -  J-'J,  '  Vs'k  ( 2575 

The  A  t,»c  la  this  case.  This  fact  effects  toe  Joottios  of  -*.  ,  the 
point  of  rotation.  If  the  section  1*  fxxsd  and  del  >rsed,  f..  ’~a  he 
arbitrarily  located  and  the  S'*  are  uaiayortant. 

When  the  section  is  sovnti*  and  is  deflected  but  not  def  need  by 
th»-  degree  oT  freedca,  determine  and  submit  t*  and  S»rk  («ttt/%ut  primes). 

The  ayrbols,  data  to  be  submitted,  computations  end  equations  used 
ij  the  SL?  are  given  tr  Appendix  II. 
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?or  sufficiently  aaalu  up liludes  of  woti —  tfc#  dynamic  effect*  or  **  i 
•io*hir.*  of  f-ci  in  a  partially  filled  tank  bare  been  analysed  la  teres  <tf 
the  r.-Atural  faude*  *cd  frequencies  of  the  easUL,  fr*s-surffc<'.»  cACilletloce  of 
tb*  fu«l.  Formulations  for  the  tolutloa  of  this  proolea  ut  pz reeat ed  for 
rectangular,  cylindrical,  and  spherical  tanka.  In  *a»v  -  cf  Appendix, 
natural  nodes  and  frequencies  are  presented  for  rec ‘.angular  *»*•  cylindrical 
tanas  la  either  a  vertical  or  a  borlxcstal  poeitlon.  la  Part  11.  ac  approxi- 
cate  procedure  Is  established  for  with  rectangular  ana  cylindrical 

tanka  that  are  neither  vertical  nor  hosts octal.  Part  111  baa  to  do  vitb  the 
treatcect  of  the  beading  node  shape  of  an  equivalent  vertical  t*.at.  Part  IV 
of  this  Appendix  concerns  the  inclusion  of  fuel  danjdag  la  the  SLP. 

PART  I 

Toe  formilae  presented  la  this  part  for  rect secular  and  cylindrical  tanka 
apply  only  to  vertical  or  horizontal  positions.  Those  presented  for  spherical 
tanka  do  not  need  to  be  qualified  aa  to  poeitlon. 

In  order  to  tUplify  the  problec  to  avpoint  where  convenient,  explicit 
solutions  could  be  obtained  in  eoet  cases,  a  surf>sr  of  assumptions  were  cade 
concerning  the  nature  of  the  fuel,  the  cotton#  of  the  fuel  and  the  shape  of 
the  tank.  The  fuel  van  attuned  to  be  noo-vltcous  and  incompressible  and  all 
tank  cottons,  except  those  normal  to  the  mm  free  surface  of  the  fuel,  vers 
restricted  to  snail  accelerations  «u*d  perturbations.  Although  the  non* viscous 
assumption  has  bean  cade,  a  damping  factor  will  be  Included  in  the  final  SET 
Structural  Loads  Program  to  accomt.  for  the  fuel  viscosity  cad  the  use  of 
baffles.  It  should  also  be  noted  that,  la  the  final  program,  provisions  are 
cade  for  auanlag  any  combination  of  rectangular,  cylindrical  or  spherical 
tanks  for  cult 1 ale  teak  vehicles. 

1.  Rectfcftcilnr  tank  -  la  the  case  of  a  rectangular  teak,  a  spring  csss 

mechanic*’  — _ ilccf  •  t _ ..  _ The  *  ^i*ace 

below  are  for  a  system  consisting  of  a  fixed  mass  K  aad  a  set  of 
undated  spring-masies  %,  so  ecnc trained  as  to  cove  only  parallel  to  the 
bottom  cf  the  taeh  cad,  in  the  case  of  the  borisoctal  rectangular  tank, 
parallel  to  the  XS-plme.  Cm  origin  of  the  axes  is  located  at  the  center 
of  gravity  of  the  undisturbed  fuel  with  the  fixed  moss  k  located  Z.  and  the 
iprlM'VMiH  at  Zt  oanetralned  by  spclags  with  stiffness  Xq  for  the  s^ 
code.  Shown  below  Is  the  horltoetal  recteegala r  tank  undergoing  longitudinal 
oscillations,  or  oscillations  la  the  2X-pl&ae.  The  notion  la  the  £6-plaa*  is 
aesuend  to  he  the  eamc  regardless  of  the  T  location.  The  following  equations 
have  been  devrloped  by  Orehea  la  Reference  (i). 
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Definitions: 

a  *  tank  length  parallel  to  X-axis 

b  m  tank  width  parallel  to  Y-axis 

c  -  tank  height  parallel  to  3-axis 

h  -  fuel  height  parallel  to  vertical  axis 

5-  -  /’abb  ~  total  fuel  aass 

S  -  fuel  aode  index  »  1,  2,  3#  —  ^  -Crdenotes  the  nuabor  of  nodes 

selected  for  use. 

f  ~  fuel  density 

2  -  acceleration  cf  tank  normal  to  ae«n  free  surface  of  Ael 

~r  «  g!U  m  total  fuel  weight 
-r  h/a  ■*  tank  aspect  ratio 

Isv  m  accent  of  inertle  about  Y-axie  if  the  fuel  were  solidified 
Ijy  m  effective  nenent  of  inertia  abewt  the  Y-nxio 
u>,  -  frequency  of  the  sWl  node  of  free  surface  oscillation 

Equations' 

U*  =[*(«' !)  £  TAHHUis-r./r^]]^ 


mJMill&zMisL. 
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7  '  n  -  nTAyHu2.S-Qg  \ 

^  ~  ~  “  fas-T)iiu 
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?--  "  vX  ^s?< 
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r  -. /  +  _Zfc&__f  t4wg.-s-,i?  •; 

Ltrl  '  M-I-,  ^  Kfi,  ►?j'rT^  4,  fiS-iTS 

hr  - 

In  the  case  of  a  horizontal  r-.ct  angular  tank  undcrgei-g  lavera^  o»cixlm.*^as 
in  tiit  YZ-plane  the  definitions  ard  equation*  are  unchanged  except  tha*  li*.* 
momenta  of  irertie  are  now  If*  and  and  the  tank  act/ect  ratio  now  hecocea 
r2  «  h/o  witU  t  tx 

£1)5  .  ($-(as-l)  -£•  TAWHiav-  >Y.T  fj.  ]]£  onJ  Is*  lir  '^P 

If  the  tank  i»  now  rotated  so  :  tho  X-axis  is  vertical  with  toe  fuel 
oscillating  parallel  to  the  XZ-.  .  me,  the  value  of  «?■/■  cbh  and  the  tank 
aspect  ratio  becomes  r,  -  h/c.  '’he  acment*  of  inertia  arc  still  taken  about 
the  y-oxis  and  the  equations  for  Zs  and  Z  ore  the  sane  except  they  become 
\8  and  X  distances  and  vk  rj  instead  of  r>.  "he  equations  for(os  and  I?y 
become* 

<ut  -  TahhU1^-'^ 

T  -  M 

iSy 

Vith  the  X-exis  vertical  but  with  th»  -3CiL^-*’ionr  parallel  to  the  XT-plane, 
*ip  t./’ebb,  the  tank  aspect  ratio  rj,  m  h/&»  tne  accents  of  inertia  ere  Ijt,  and 
I3J  and  the  equations  for  ft?  B  ru*xi  I8-  beccoe: 

tos=  [?-<is-0^  TAvnUi5-i)rrrf])l'z 
r  -  M-  itii.1 

*5*  n‘r  It 

Ibis  ccrr>letci  the  specification  of  the  equations  for  the  spring -mass  analogy 
for  rectangular  tanks  in  a  horizontal  or  vertical  orientation,  7he**efor'*,) 
the  angle  the  X-axis  uakes  with  the  horizontal  determines  whies  set  of 
equations  wore  accurately  approximates  the  situation. 

?.  Cylindrical  Tank  -  The  formations  for  the  cylindrical  tank  are  not 
nearly  as  straight  forward  at  vert  those  for  the  rectangular  tanks.  Three 
different  jetbod3  ha/c  been  used  to  define  the  fuel  action  for  the  different 


tarJ  orientatic-s  As  extensile  literatus  survey  indicated  *aav  Tor  the  case 
of  -»s  Illations  in  a  Uoi  itootoA  cyMn^ri -si  ■* ».  V  at  '<.”*■*  approach 

ccsild  be  use  I  as  in  the  case  *f  rectanciar  tanks  Although  no  development  of 
the  equa.vic*.c  cculd  be  found,  "e ferer.ee  (2)  suggested  tbt  tie  natural  fre¬ 
quencies  fo.  tic  norizcntal  cylindrical  tank  are: 


to;  TMH  Zfh]  ** 


•where  s  -  1,  ?,  3,  Usr 


Scsp&rlrg  thii  equation  th  the  corresponding  i. ecuency  equeMor  for 
rectangular  tanks  indie/, ten  that  the  cylindrical  tank  aspect  ’--tic  is  r  -  h/l. 
Making  like  cocparirons  tX.c  follow i*v>  development  is  '  ugge stod 


<z 


i 

9 

E— - /  — 

definitions: 

i  »  tank  lrngth  parallel  to  X-axic 
R  •  tank  re ilus 

h  «*  fuel  height  parallel  to  vertical  axis 
Jl.  •  total  fuel  sass 
»  -  1,  i,  3,  —  tw. 

r  «  h/l  •  tank  aspect  ratio 

Squaticnn: 

l-fpTAVHUm-)]**- 


~6l~ 


-  ^  -«5»w£*4j3B 


-  i 


-^r5>  TTv- 


TT  2>s*s 

5*1 


^■FY'UyJ1'  ttf»- 
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rhe  second  method  to  be  used  on  the  cylindrical  tank  folio vs  the  formu¬ 
lations  of  J.  '•/.  Miles  found  ir«  Reference  (3)  for  an  *t. eulnr  cylinder. 

In  thl3  analysis  the  potential  end  kinetic  energy  expressions  are  derived 
with  allowances  mede  for  tan!:  flexibility.  First  the  potential  energy  (u) 
and  the  kinetic  energy'  (T)  expressions  are  stated  and  then  the  potential 
energy  coefficients  and  the  inertia  coefficients  (u^j )  are  defined. 


Definitions: 

i  -  j  .  a  -  1,  2,  3,  —  «r 
q  (t)  m  Generalised  coordinates 

ql(t)  —  a  translation  along  e  «■  o  q 

q2(t  j  —  a  rotation  about  the  ccntroidal  axis  f>  m  2 
q-(t)f  (X)  *  a  single  bending  displacement  along  6  m  o 
q^S^g+^r,©)  «■  sloshing  displacements 
f (X)  m  bending  mode  shape  of  tank 
f(X)  -  df(X)/dx 

l^a+jfr,©)  m  sttl  mode  shape  of  fuel 
M  m  2f TP  R^b  m  total  mase  of  fuel 
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Definitions  ( "'cntinued) 

s  «  index  Indicating  fuel  slosh  tcdcs 
l>  *  one  half  fluid  height 
.1  ■  tank  radius 

a  *  *cotleratioa  of  tank  along  X-axis 

'  stil  *cro  of  the  !irst  derivative  of  the  FLnetior  of 

the  first  order  and  the  first  kiwi. 

(ySx  •  1.0!>li9,^2  ■  .  8.53631, ck  *  11.796001 


•  1  If  1  -  J 
0  if  1  i  j 


liquations 


Peference  (3)  defines  the  potential  energy  coefficients,  as  shovrt  helov: 


4t ,,  ~  'K,t  ~  '&U  =  '£/$  -  ?  s 

*  *  r 

S‘xf'w“y 

^*♦•>4*1  “  it  (^1  ~0 


T-t^/Wg^y 
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rt’,j  "j*,  are  4c*vx£  cv 


^.,  -  Aj 


05}  • 


>‘,1 


!'*'-1 '-  *"j;  r  M j  rc  -  J—  tff-)  -  7  (-cj)  i 
^/5T3j(  = 

™  -  Mfb1-.  3Z'\  ^SMR3  “  TAf.-n(^ 

T  -r^-s-  l  -BTTepfr 

-  *%.  =  +*{ii  xfuUi 

+[?(»+?(-*][$-$**£  I  ^ 

*>*,*»  -  =  0«  [ ian*eJM-,] 

■™3} ~ Mr'  ->■■££'% Fs*~  nr  I 

+  ^ml[?U) -fh]  T^r-Z  L^Zs.[j-)i"t'(b)rnv]r3 


P  £i 


/ 

^3  -  ^3-ft  V*- %1  [CSCH  (I^Jj^MWH^i)] 
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1 2  2ji<3  Io  are  certain  Bessel  functions  c»r  c<»rt.'ii  for  -a  -i  functions 
of  com?  related  kind. 

C0S^'r'irs)^x  whsn  s>°  • 

k  xtt  lt 

These  equations  have  been  used  in  the  experimental  analysis  of  Reference  (-). 

In  this  report,  resonant  bending  frequencies  and  mode  shapes  were  d'^reined 
experimentally  and  were  shovn  to  be  genera’ lv  in  cgrecnent  viph  the  ->rcticol 
predictions.  The  deferences  vre  attributed  p.-.marily  te  variations  of 
actual  wfie  shapes  from  those  *ssu.med  in  the  theory. 

"he  third  method  to  be  used  on  cylindrical  tanks  -ras  formulated  by 
3.  ’iriiar.sky  xn  Reference  (5)*  This  is  elso  the  method  to  be  u#ed  ior 
spherical  tanks.  In  this  report  a.,  ir.tecraiequatlon  approach  is  u3cd  and  the 
method  of  solution  developed  for  the  fi»*st  three  fuel  slosh  node  ,  «uich 
as  indicated  in  the  literature  is  a  sufficient  nuiber  of  nodes  for  cost 
rract Leal  problems.  The  tank  orientation  under  ccrr-Ideration  is  c  horizontal 
c>li..dr*cnl  tank  undergoing  lrtej-al  oscillations.  In  this  case  the  ^.urali-cJ 
zoor< i-atc  denot' nc  ration  slor._.  the  Y-axis  is  q1  and  vain  the  fuel  sloahitg 
gencnlixcd  coordinates  are  ^**3.  dimensions  of  the  q6*’  in  this 
dove1  spnent,  hove  vs  r,  aix-  (length)"  rateer  then  length  a3  ir.  t  \e  case  of  the 
pro vi ran  fC'cralircd  coordinates.  '~ic  slosh  heicht,T  J+3,  at  the  side  cf  the 
tank  can  be ^exp cessed  ns  a  function  cf  the  by  the  folio  .Inc  relation, 

m  q^^Clj  -a+j/f*  It  should  be  noted  that  in  this  a-iol/axs  tank  beudi.g, 
ic  ignored  anu  that  with  the  .icn-viscous  assumption  rotation  of  this 
cylindrical  tnrJ:  tmi  rotatic  i  of  th*  "pb^rjr*’  tnrfc  •'eed  not  b*»  considered. 
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Definitions : 


q1  m  >i  ‘••vctiiolitio'j  tlon;:  the  a -axis 

f  tame  radii’ 3 

c  =»  fuel  Vieichc  purs-.stov  =»  -1.0  for  empty  tarJ* 

»  sin  jQ  ■=■  0  for  half  full  taul* 

a  +1.0  for  full  tank 

&  ■»  si-.r^o 

tt  m  COS  0 

jNs+3  =  frequency  parameter  «  ^  Cl)  s+3 

s  «s  fuel  slosh  node  index  =*  1,  2,  3;  —  Ur 


i  -  tank  length 
*r  «  fuel  density 


tanl  acceleration  alon;,  1  -axis 
total  mass  of  fuel 


equations : 

V-rtlUt 

O  S*l 

U  =  AnAf')' 

K  $,/ 


T-kMf(jo*+ l  <"w  *•»&*)' + <*u.r 
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Z\k  rv:-Uirs?nsionaJ  30.  >1  pareceters  and  \+-$  aloug  v;  +  h  jXs+3 

;t  pr^cc  ved  In  -igures  (k)>  (5^  and  (6)  respectively  for  the  firs*,  three 
aljj>h  codes  furc* Ijne  of  the  ftcl  height  par  octets*  e.  It  should  be 
rv  ted  that  for  values  '  ♦*  e  *  .H  to  e  •  1.0,  the  curves  in  Fhijures  \6)  a*xl 
(9)  tend  to  infinity  as  thr.r  approach  e  ■  1.0.  To  id  ,#,io  type  of 
solution  near  4  *  1.0,  tt.e  curves  have  beer  tsci  v*  Intfre...  t  e  *.  1.0  to 
provide  for  an  approximate  out  finite  solution  fer  the  full  wnic. 
this  reason  the  solution  of  the  equations  for  the  r«’e;ly  fuj.1  to  full 
cylindrical  and  spherical  tanks  oust  be  used  with  cautior. 

3.  Suh^rlcal  Tanks  -  ?o l loving  the  ease  method  anrf  definitions  as 
used  above,  the  solutions  for  the  spherical  tank  nay  be  obtained. 

Equations  U-  (“■£"*'  £  CO4.,  C„,  (g"*”) 

»  »-f 

T=  WS» 

Oi*,.,  D«.  j*+’ 

T  »  0*  «•  *  L  [VCO  v_ „ , 

As  before,  the  values  of  the  noadiaerslonal  sods'  p^aoeters  0*43  and  0,«3 
along  with 4!\7+3  aro  plotted  versus  e  in  Figures  ,'t),  (8).  end  (9; 
respectively.  As  discussed  previously,  the  solutions  of  the  eolations  -"e 
only  approximate  solutions  vs  the  full  condition  is  approached. 
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Figure  4.  Variation  of  A  with  Fuel  Height  Parameter,  » 


Cylindrical  Tank 
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I  \  Figure  5 


$• 


..p  ..4  _a  6  .2  .4  .*  .*  1  c 

Variation  of  3  with  ?uel  H-aigSf  Po  'anatir,  «f  Cylindrical  Tank 
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«i.a  -i  -A  -.J  0  .2  -4  *  .*  1.0 

ftgur*  7.  Variation  o?Cj+3  with  Fu*r  Haight  Parameter,  •,  Spherical  Tank 
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F:^ur»  8,  Variation  of  with  Fvol  Hoight  Paraw.otor,  •,  . 

Sohorica!  Tank 


Figure  9.  Variation  of  Spherical  Frequency  Parameter 
wirh  Fuel  Height  Parameter 
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"oe  cccwtations  cabled  for  in  *M«  part  axe  '  *rr led  cut  1  -  tie  -'-Uxcie 
Pte'«l'*'l  stf  *«  ^bpresras  (vpCSS),  •Tact  runs  »- eh  the  basic  SEP 

prog --vi  -ether  then  tbs  SLP. 


Tbr  free  surface  of  the  fuel  la  parallel  to  the  horjr*  n  cujy  when  the 
tW*  tas  no  lateral  ox  longitudinal  acceleration.  it  ic  sovldpctel,  however, 
that  such  vil.  act  usually  u!  the  case.  The  »y:n*tption  no<f  being  node 
concerning  the  fuel  orientation  is  that  the  ?r*c  r.vf ace  always 

perpendicular  to  the  "resultant  tank  acceleration,"  detlued  •  the  actual 
acceleration  at  the  tank  center  due  to  the  gross  notice  of  tb'  vehicle  ainue 
the  force  per  unit  case  due  to  gravity.  The  pertinent  angle:  ter  the  task 
orientation,  therefore,  are  not  angles  usually  defined  as  the  tank  o- 
vehi.de  pitch,  roll,  and  yuv  armies;  they  are  the  angles  between  the  body 
axis  system  and  the  resultant  acceleration.  This  aesas  that  at  any  Instant 
of  time  tb*  resultant  tank  acceleration  crust  first  be  found  and  then  the 
free  surface  of  the  fuel  set  perpendicular  to  it.  Since  the  fuel  slosh 
equations  presented  in  Part  I  are  valid  i%r  only  vertical  or  horizort-al  anas, 
the  tank  vails  must  be  set  perpendicular  aar  parallel  to  the  free  surface. 

As  this  is  dene,  the  real  tax*  dimensions  in  the  body  axis  systee  are  replaced 
by  thoie  of  a  different  "but  "equivalent*  tank  of  the  some  volume.  This 
equivalent  tank  is,  therefore,  a  tank  vhose  dimensions  and  orientation  are 
a  function  cf  the  angles  the  real  tank  makes  with  the  resultant  tank 
acceleration .  As  the  real  tank  for  example  pitches  from  0°  to  90°  f  the 
equivalent  tank  concept  provides  a  continuous  transition  to  classify  the 
tank  as  being  either  vertical  or  horizontal.  The  tank  geometrical  corner 
was  chosen  as  being  comnon  to  both  the  real  tank  and  the  equivalent  tank. 

Tfcs  equivalent  tanx  concept  is  by  no  -cans  as  exact  representation  tut 
does  give  an  approximation  of  the  real  situation.  One  very  significant 
parameter  to  fuel  sloshing  is  the  length  of  the  free  surface.  The  equivalent 
tank  concept  permit*  the  free  surface  length  to  increase  or  decref^e  as  it 
does  in  the  real  situation,  but  only  apprw'^^tes  the  actual  freu  surface 
length.  This  concept  also  simplifies  the  computation  of  the  moments  and 
products  of  Inertia  and  the  C.C.  of  the  fuel,  as  the  fuel  changes  its  gross 
position  in  the  tac*-  ’ue  tc  the  gross  motion  of  the  vehicles. 

Consider  now  the  problem  of  obtaining  the  equivalent  rectangular  tank 
dimensions  and  then  the  moments  cf  inertia  ~nd  C.C.  ccf  the  fuel  in  the  equivalent 
tank  as  If  tSe  fuel  were  solid  fled.  As  sb.  m  In  Figure  10,  X  is  the  unit 
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»cceler.ition  vector  of  the  tank.  One  corner  of  the  real  tank  la  choeen  a c  the 
origin  of  n  right  handed  triad  having  unit  vectors -yi,,  i* ^  and'/7J,  pointing  along 
adjacent  edges  of  the  tank  and  chosen  no  tliat  the  angle  between  i3  and  or,  is  not 
leas  than  the  rvngle  between  Ay  and  either  ->ii  or  .  Jbe  «'  urtex  for  T.nc  comnon 
oriel**  of  these  vectors  is  chosen  so  that,  these  angl.-s  arc  not  greater  than  rr/?- 
The  uni  t.  >/;i’tor  &  i r.  defined  as  a  unit  vector  perja  ndicul:*;  to  -W[t  lying  in  the 
p.l*>\ .  zi  . , j  and  •>»,  ,  and  muking  an  acute  angle  vith.-fj  .  The  vectors  J£,  ,  2  t 
and  2J  txra  the  vectors  defining  the  real  ton!.  .  Lze  and  orientation.  TV  i^icgn.  tu* -es 
of  these-  vectors  (not  necessarily  resper-ti-ely)  JL\  ,  .4,  and  ,4jarc  the _l-*ngth,:  ui 
tne  sides  of  the  tank  in  the  direction  of  the  unit  vectors  -h,  ,  -pftand''*H  respec¬ 
tively,  as  shown  in  Figure  10.  The  unit  vector  raid  the  angles  0(  ,  and  Q%  can 
be  defined  as: 


0t  *  ARC  COS  (/3-e) 


~  ARC  Cos  i/7-(Fyfn, )* 

'! 


9,  -  ARC  Cos  (fjy-e) 


=  ARC  Cos  IJysrTjSSceJ 

where  Q,  and  0^  oust  be  pos-tive  acute  angles. 

The  dimensions  of  the  equivalent  recteu-jolar  tank  can  now  be  obtained. 
Referring  to  Figure  10,  the  equivalent  tank  may  be  thought  of  aa  the  tank 
obtained  by  taking  the  real  tank,  with  7$  coincident  with  one  of  its  edges, 
and  then  adjusting  the  real  tank  dimensions  to  the  equivalent  tank  dimensions 
as  the  tank  is  rotated  first  through  &,  ,  and  then  through  .  This  then 
replaces  the  real  tank,  which  is  actually  in  the  position  described  by  &, 
and  £kbut  with  its  sides  not  perpendicular  to  the  free  surface,  by  an  equiva¬ 
lent  tank  with  the  same  volume  and-  approximately  the  same  free  surface  length 
with  its  sides  perpendicular  to  the  free  surface.  Defined  below  are  the  equiva¬ 
lent  t-«»rfe  dimensions  C, ,  C« ,  and  Cy  in  terms  of  the  real  tank  dimensions  4, , 
4^and  t!a .  The  intermediate  tank  dimension  C  is  defined  as  the  length  of  Cj 
after  the  tank  has  been  rotated  through  Q,  but  not  through 
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c-  JijTl  Tm[o,  r(i-  arc  Tan  -/Zyi't] 

C,  -  ~J'A c  CoT[el+(i -  aactan  vc/Zr  i 

-  '~,7^T,  cor[9,+(i  -  *&)  arc  Tan  vCT^yj:  ) 

e,  -  -JTTc'  Tan  la,.  ^a-i±-)  pK  -tom  vT/ts) 


The  unit  vectcrs  giving  *n‘*  » **w  —  •  ■  ■  ->  ore. 


Equation*  fer  the  f-?l  counts  of  inertia,  J,  ,  and  J, ,  *s_lf  the  fuel  were 
solidified,  token  obouv  the  fuel  C.C. ,  and  the  C.C.  location  Zp  of  the  fuel, 
oeasured  froe  the  ecuivolcnt  tank  center  alone' the  A  axis  are  sh(  m  below. 

The  total  =ass  of  fuel  in  t»*?  tank  is  and  the  fuef  height  along  ^  i6  K  . 

T,  = 


I,  =  -k  Mf[(c.)l+(co‘] 
if  ~  _4r  ■< } 

K6  M  f  /  fo.C^ 


An  approach  dollar  to  that  used  for  the  rectangular  tans  is  presented  for 
the  cylindrical  tank.  There  are  two  ca^or  differences  between  the  equivalent 
rectangular  and  cylindrical  t._Aks.  The  cross  section  of  the  equivalent  tana, 
taken  perpendicular  to  the  resultant -acceleration,  is  always  rectangular  for 
the  equivalent  rectangular  tens.  This  cross  section  for  the  equivalent  ~ylln* 
drical  tnnk  any  be  rectangular  cr  circular  demanding  on  the  ar^le  betwew  the 
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res-itm-i  acceleration  and  the  rear.  'onk  length  "jetor  If  *h*s  '•oss  section 
is  circular  the  equi/alcnt  tan*  is  considered  as  being  vertical,  and  if  rec¬ 
tangular,  the  eq*ivalent  tans  is  considered  as  being  horizontal  tic  second 
dlfferet  n«  betvejn  the  equivalent  rectangular  and  cylindrical  tt../.  con-crns 
tr«.  an-*l»s  act. d  to  specify  their  orientation.  For  the  equivalent  rcv.-u~tg-t.-ar 
t— 1»  two  «n«Ucs  were  needed.  For  t hr  equivalent  cv._nd:ical  tame,  the  angle 
ueOiC*  ate  resultant  acceleration  and  the  tank  length  vtr*or  j£  is  tne  only 
i  *lea  to  specify  the  tank  orientation. 

A.*  eho-n  in  Figure  11,  t  is  the  cylindrical  tank  length  vector,  is  the 
reel J tent  acceiiretj.cn  unit  vector  ar.d O  U  *he  angle  oet  ••■».*»  then.  Defining 
<M  the  tank  length  uni*  \ actor  then 

I  *  I  A 

y  »  -f3  ’A 

9  *  fitCCDS  lUl 

If  1^1  it  greater  than  equivalent  tank  1b  vertical.  Defining  -£cnd  "P, 

ns  tne  real  tank  length  and  r-oius  respectively,  and  ^and  H¥  as  the  equiva¬ 
lent  tank  length  and  radius  respectively,  then  the  relation  between  then  in: 

X-y  -  Mzlh  T/°v{e+(/-  $£)  T?>»  vCt/ln] 

?y  =  VvxZu 

h  =  Mr/rrf-Kf 

Equations  for  the  cooents  of  inertia  of  the  fuel,  a*  if  the  f-el  -ere  solid*  - 
fied,  taken  nbout  the  fuel  C.G. ,  J,  ,  t  **  .  and  the  C.O.  location 

of  the  fuel  oeaaured  along  ^  froa  ♦be  teak  center  are  shown  below. 


■ij  = 

Wf  -WM) 

If  ({**>  1-0,  the  components  of  the  unit  vector*  %  end  giving  the  nev  cilrec- 
tioac  are: 


.  .  M; 

r*  ’  o  A\ 

y'  , 

z  ■f'l-aiy1 

•JT-(tl)- 

„  -.u,  t  .e  components  of 

^  and  -7^  are 

if 

1 

1 

A‘  X\-W'  ■ 

A>  u£ 

sl  .  -ti 

Jl\  --  -?.}  f 

V  t-  mr 

-  JiA'-Jli1 

V  /-(u)*- 

V  |-W« 

If  1*  less  than  or  ecc-ai  to  the  equivalent  cyAindri-al  tent  is 
horizontal.  Defining  X^-id  7f*  -*»  the  equivaxenv  ronzonzel  teni  lengtn  and 
radius  respectively,  and  expressing  thtc  as  functions  of  the  real  t*ck  di sua¬ 
sions  !■  and  7?  gives 

•£»  -  VeTW  coT[e  +  (i-p)  wtrwi^] 
fix  -  V/PV-/'4-* 

In  order  to  obtain  the  expressions  Cor  the  cooents  of  inertia  and  C.3,  of  the 
horizontal  tea*,  the  angle  $  ,  shova  in  figure  \S ,  aust  be  found.  The  equa¬ 
tion  relating^  to  the  fuel  enss/f^  is 

Mr  3  ftx  (  %  +  B  +  *IIJB  CcsS) 

[(Ml/f-i-x  K ")-  Jgr]  -  B  +  b'~0CO!>0 

Xcv  let  Nevtoa's  method  can  then  be  used  to  find  B_ 

Defining  f(fi)  end  shovr  belov: 

f(fi)  *  B  +  SIH0  coi  £  —  CH 

f '(.»)-  e  cos'B 

Let  B  uf  U.e  inltltl  estimate  of  $  aod  calculate^ 

A*  £<h 
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V*»  sx 


'■'-&*-  %$)■ 

The  quf...  titles  J  f^j!  “MM)  l  can  then  he  tested  to  determine  if  they 

are  both  lea*  than  (l  x  10*7 ) }  the  arbitrarily  chosen  degree  of  accuracy.  If 
this  is  true  then  /6  m  B,  •  If  the  desired  degree  of  accuracy  baa  not  been 
obtained,  /?  la  used  aa  the  next  estiaatc  of  4" ,  and  a  ,£^  ?.iet  he  calculated. 

The  teat  la  made  again  to  determine  if  the  desired  value  o tB  nos  been  obtained, 
end  if  not,  the  iteration  ouct  be  continued  until  the  deal  red  conditions  are 
uetisfiad.  Tha  fuel  moments  of  inertia,  C.G.  and  height  can  '■•“'.an  be  obtained 
as  functions  of  /$• 

2F  =  -  ZSnSJ&HB _ _ 

F  3[‘%+£+ SIKH  Cos  BJ 

«T|S  +  ZF(?1tH  S\HB~2F)] 

!T«r 

J» 1  +"{t  si»s] 

2,  =  ^ 


Tha  spherical  tank  dimensions  do  not  neea  adjustment  because  for  any  tank 
orientation,  tha  free  surface  length  will  remain  unchanged.  The  orientation 
of  tha  free  surface  within  the  tank  will,  however,  change  positions  In  the 
tank.  The  angle £  for  the  spherical  tank  is  defined  in  figure  12,  and  the 
came  iteration  method  as  described  previously  must  be  uaad  to  solve  for  A. 

The  following  equations  must  be  used  for  this  Iteration.  ~ 

Mr*  fr-fS’f^vsiwa-^.  &!#*£) 

I  Trfc r*-  H  * 
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{(#)  =  Siv&-  ^  S;ij]s  -  C 


f*(8)  -  Co5  %& 


P.*C 


Using  the  value  of  &  obtained  fro*  the  iteration,  the  following  equation! 
for  the  accent  of  inertia,  fuel  height,  and  C.C.  can  be  solved. 


2 

~r 


-  3P  Cos'S _ 

dfs,/  S/V*C) 


f»  *  ft  U+s/v/t) 


T,*  rt-  frMrfe'Vt+fZr  n$'U/e-SK) 


J,=  ftorW'-irffstv*) 

3he  coepcoectt  of  ^  and  ^.giving  the  new  directions  are: 


I 


FAirr  m 


Tnnf!C>  as  the  length  of  an  •equivalent*  vertical  cylindrical  tank  Is 
sort  likely  to  he  different  fro*  that  of  the  real  tank,  .<ad  since  the  tending 
node  shape  7  U.  is  given  for  the  length  of  the  real  t->r-ic  but  mist  be  applied 
along  the  length  of  the  .equivalent  tank.  It  le  necessary  tc  find  soae  vey 
of  ad«pti«g  the  use  of  7  l«l  to  a  changi*'*  tank  lezgVi. 

* 

Vo  a-.^aplieh  Ibis,  a  jpev  arpae at  % «  Is  Introduced  which  dees  not  very 
with  the  tank  length,  and  r  is  given  as  /  i*  i ; .  Separate  emsUer^  b.-.  cunt 
he  given  to  the  two  cases  ui  *  */VF  and  u;  <  -  'A/f,  where  e»j  equals 
CCj  •;  and  for  a  vertical  'equivalent'  tank  Suii>  : /VT  •  *'i  i»  defined 
as  distance  along  the  axis  measured  froa  the  center  of  cylindrical  tank  >  , 

nondlnenslonallsed  with  respect  to  the  tank  length..  Xa  use  /1th  the  real  tank. 

X  •  I-  positive  In  the  direction  so  chosen  la  connection  wit"  the  sahnlssion 
of  as  .a  to  the  VFCS2.  In  use  with  a  vertical  equivalent  tank.  x'i  is 
positive  la  the  direction  at  the  'resultant  acceleration"  If  tao  >/Ve‘  ,  and 
positive  In  the  direction  opposite  to  the  "resultant  acceleration ’  if  u.;<-  -»/vr. 


Far  the  suhalsslcn  of  data,  we  note  that 

As  i  is  the  length  of  the  real  tank,  that  Jr;  X  >  is  actual  distance  ■ensured 
in  feet  fro*  the  center  of  the  real  task,  that  -  <  a  4  ,  and  that 

X  ’  M>  r  is  #.  A*  * 

Ve  next  consider  the  vetheaetlcal  relations  coons cted  with  the  equivalent  tack. 
When  U  j  >  I/V2T  ,  these  are  as  shown  below: 


X ;  *  lrv,  X  ;  t  (,  Irwj  ~hi)  /<L 


& 

hi 


1 1 


Ax i  -  lrvt  i  X  i 


X  ■  *  X;  -  C ir v; - V;-- )  U 
kvi 

'fc,  (a  '  ft)  -  s 

Whew  X;  «-  lOi  ,  *'{  «-X 

9  * 
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When  Xi 


*  S,  i  < 


A  f 

•vi  * 

^  L 


*  j-i 


J2j  _  1 

e  J?  r  /  i 

•'i  1.0s  |  =>  »  -t—  v  a 
L  L  «  i 


v  ''  '•  1  r. 

■<  ;j  «  Ai 


£  •  f/(-i) 

t  -fUt-i) 


&1  '  *"/.? 


wl5*n  u  i  "c  ~  »  tho  pirtir-snt  relations  are  as  foil  own: 

X  i  *  l irvj -Vi.)/*  -  irw}  X; 


iX;>  ~  bvi  &*\  j  (\\\\\\\N^  T 

iii  ,  Jt?  (4-x^*i  !  K  ^  '  \\\\\  I 

«■  Tfc  1.1/*  X*  '  pjj-  H-v.-W. 

». v\\VA\\  I 

»  S  Wj  /I  V;  V  „■  ,  wj  * 

«  "  t.i  U  *;/ 

WW*vt  X  i  *  -  ,  X ;  *  g  ^ _ 

wk.  1 

'  Jt.i  f  *  *“£  _  Jbs  f  ^  fix'. 

F«  -  -  “4  MXl  *•  4-^t  ' 


"Equivalent  Tank' 


pwi?  r/ 


Aa  indicated  In  ec’iatioiis  (§7)  and  (95),  provisions  L/av?  been  itaae  to 
incJjide  in  tbe  SI?,  iurov  ledge  of  whet  numerical  values  to  use  for 

''■>  6  j  )  la  iaportant  to  a  careful  investigetlce.  of  fuel  sloshing  In 
t>  vehicle  *3  * light.  Aa  extensive  literature  search  found  that  very  little 
fuel  <us 3*1.^,  data  exists  except  for  upright  cylindrical  tanks.  Reference  (6) 
did  present  tbs  eqiatton  below,  which  can  he  u.~a  to  obtain  the  logarithm!- 
decrement  b ;  for  an  upright  cylindrical  tank  as  a  function  of  the  kineuttic 
viscosity  v  >  the  fuel  height h  ,  the  uccejerstloc  due  to  gravity  g,  and  the 
tark  radiu*  R: 


6  *  C»  ♦  g  ( « -  h  /a)  cich  is  .<>%»>  /r  ^3 

t»<*W  l5.*4h 

*fcis  aquation  is  for  a  tank  with  no  taffies.  *fc*rt  of  the  other  references 
food  were  for  upright  cylindrical  tanka  with  various  baffllrq  configurations. 


of  the  scarcity  of  data  cn  fuel  Anqing,  no  equations  such  as 
the  am  Just  given  (which  is  at  Halted  applicability)  ere  Mqlcyed  in  the 
81?.  lather,  it  Is  left  to  the  user  to  detsratos  la  hie  awn  way  constant 
values  at*')  tar  sobelttal  aa  input  to  the  grqgrsq.  Aa  lcag  aa  the  schattted 
values  at  <f")  are  greater  than  sere,  they  will  at  laaat  prevent  the  infinite 
continuation  at  whatever  fuel  slosh  nodes  are  excited  ty  the  action  of  tbe 
vehicle. 


Symbol  o 
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SYMBOLS,  DATA  TO  BE  SUBGTTED, 
rr£zzTPL’?:rxz  uKn  c^iKiluio  USSU 
IB  THE  STRUCTURAL  LOADS  PROGRAM 


components  in  the  y  coordinate  syutoe  of  the  linear  eccelero  r,ion 
o'-  the  vehicle  at  the  origin  of  the  vehicle  axes. 

SIAR7T 


£  components  in  the  y  coordinate  system  of  the  linear  acceleration 

E*J»  of  the  h-th  particle  of  Lhe  i-th  section  due  to  era.  tic  deformation. 
APAEHT 


Ar  components  in  the  y  coordinate  system  of  the  linear  acceleration 

Riil  of  the  h-th  particle  of  the  i-th  section  due  to  rigid  motion. 
ADAR3F 


Ara 

;  static  aerodynamic  terms. 

M 


S(-  nondimeurional  fuel  slosh  modal  parameters  for  lateral  motion  of 
horizontal  cylindrical  tanks. 


A  a  quantity  used  with  vertical  cylindrical  tanks. 

TAAPS  _ 


sectional  aerodynamic  shear  force  terms  for  rigid  vehicle,  referred 
*  to  vehicle  axes. 

SAAPTT 


sectioral  aerodynamic  bending  moment  terms  for  rigid  vehicle, 
1  referred  to  vehicle  axes. 

SAAPTr 


A  ci*,  *  A 


components  in  the  y  coordinate  system 
of  the  linear  acceleration  of  the  h-th 
particle  cf  the  i-th  section. 

ADAXHT 


*  ■ 


/ 


inertia  terns. 
S FLAPS 


c:--\o to  iu  the  V i  ixr.ate  system  .t‘  the  given  mode  of 
••ibretion  in  degree  of  freedom  fc  before  balancing. 


cosiponents  in  the  system  of  the  i/’ven  nr-tial  'incur  velocity 

with  respect  to  q*  of  the  ooint  of  rotation  oT  .  o.-sbie  section  i 
reiative  to  the  vihicie  before  balancing. 


aA  nodal  functions  of  cylindrical  tank  aspect  ratios. 

TAAH 


CUluti  nodal  functions  of  rectangular  tank  aspect  ratios. 

1AAM 


_ .  nondisaensionel  fuel  slosh  nodal  parameters  for  lateral  motion  of 
horizontal  cylindrical  tanka.. 


ft  aerodynamic  stiffness  terns. 

j  K 


_  /rib 

D  sectional  aerodynamic  shear  force  terms >  referred  to  vehicle  exet. 

*'  SABPIT 


»  rwt 

0»  sectional  aerodynamic  bending  aooeut  toms,  referred  to  vehicle  axes 

SABFPT 


length*  of  equivalent  horizontal  cylindrical  tanks.  Same  as  in 
VPCS.  "  * 

TABH7T 


\Ls—  +rt 


<■.  ■  •* 


"***  **"  ■W-f'  CfS  W  »  *>»*  JW"  «'»■'» •»"i‘  '"«•*  '"“V 


Sane  as  in 


lengths  of  equivalent  vertl.:*!  cylindrical  tanks. 

7>C->. 

TAi>m 

r't  dynamically  balancing  rotati  rata  with  ••c.jpeet 
tof^cf  tie  vehicle  relative  to  the  vehicle  axe.'  • 

■itSBFTi’  -  sraucrjiE 

hABFTT  -  TAl.KS 


Q&>  nenditensionej.  foal  sloch  acdal  parbLX-l-.n ^  lor  i: ~ sral  tnotion  of 

spherics1,  tanks. 


cercdycamic  damping  terms. 


c 


'  rs 


tu 


sectional  aerodynamic  shear  forte  terns  referred  to  vehicle  exes. 
SACPIT 


i 


sectional  aerodynamic  bending  moment  terns  referred  to  vehicle  exes. 
SACFFT 


C  •  )  C*%  lengths  of  the  "horizontal"  edges  of  the  (equivalent) 

14  rectangular  tanka.  Sane  as  ia  VBCb. 

TACaTT 


components  of  the  dynamically  balancing  translation  rate  vith 
respect  to  <^*of  the  vehicle  relative  to  the  vehicle  axes. 
SKCTJCT  -  SSRK3USB 
TAtTCT  -  TASKS 


0^.  nor.dlscnslonsl  fuel  slosh  nodal  parameters  for  lateral  notion  of 
1  spherical  tanks. 


5 


* 


ttS*. 


■s%  ;  -s-r  3%- »  -xn  'v'.# 


n 


"r-entJ  onti  r.-.gtti f  or.  cut  to  of  inertia  of  ;tructurK  sr..* 

?■&  '‘ucl  -ihoat  axes  thru  the  venlcle  center  of  n»  3  and  par»i3ei  to  the 
/•<Mde  exes. 

'.1.77S 


O’.JK 


inflect  loro  of  the  il  -th  particle  of  the  t  -t.t  oe-tion 
•.  lr-tic  o  -fonatior; . 

AlnlTHT 


UUC  CU 


mcvJa i  product'*,  of  Inertia  wf  i.art  j.%  i*vj  vehicle . 

rvc,  V 

r  CDDK3 


d  •  dynar  ic 

J*  ?nwK 


’-ini.cir.g  "era. 


El  nui.be r  of  thru-t  vectoring  nozzle",  (or  "engines '*}. 


K03HG 


£,  rood  inertin  ts/r-. 


3r2K7S 


6  base  of  natural  -v  .'-."tea. 


6  .  components  in  the  J r  oyjtcn  of  the  vectori . 

TAESRT 


Saae  as  in  /?C3. 


r  pjj  raoaents  and  negatives  of  products  of  inertia  of  fuel  about  vehicle 

8X83. 

TAFR?S 


F  .  certain  integral-  connected  clth  fuel  aloah  ir.  vertical  cylindrical 
tanks. 


TJUfiPS 
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1  U ./* 


\M 


VV. 


H  nodal  ur. balances . 

*  ‘  br-fTlC 


h. ,  weight  of  fuel  in  "equivalent"  tank  i  •  3  :se  as  is  VPCS. 


components  of  the  partial  lines/  velocity  with  n.  ,r-ct  to  q J  of 
the  center  of  aass  of  section  1  relative  to  the  vehicle  axe..  --  value, 
obtained  after  cynaole  balancing. 

TAHKI 

SFHKI 


H  F  codal  unbalances  for  fuel. 

w  K  T  AH  FT’S 


H  F  inertia  couclir.g  terns  for  fuel. 

u  jK  *  ° 

TAHFJS 


HS-.  v  inertia  coupling  terms  for  structure. 


tr.iK.  ‘  - 

oarjsoo 


HS  nodal  unbalances  for  structure. 
*  SSH3TS 


K  ..  noneuts  and  negatives  of  products  of  Inertia  for  part  of  a  section. 


SIHHIS 


nooents  and  negatives  of  products  c- 
fuel  about  vehicle  axes . 

CFIRSS 


*  Inertia  Of  StruOtUl'c  ttiiu 


noaentc  of  inertia  of  fuel  in  ''equivalent*  tanks  about  axes  parollel 
to  vehicle  axes. 


~jT  j.,  ftoKer.ts  and  negatives  of  products  of  inertia  for  pert  of  the  vehicle. 
SIIIBS 


t  l  ■ 

®  •.  ..  t 


T 


»?•  £"lvrji^.a 

>  +  . 


accents  of  inertia  -f  sootier.  1  about  . ecu lor. a  axe.  Same  us 

\n?'. s . 

TAJRVS 


aor.cn  ts  of  inertia  of  fuel  os  if  it  vere  ocJiu  •  ■'equivalent " 

incut  tank  axes.  Sane  as  in  VPCS. 

iprrjs 


effective  moments  of  inertia  of  fuel  a  oou ;  «  ■». 

TAJFPS 

* 

components  of  the  partial  linear  velocity  uith  respe.r  uc  q  of 
the  center  of  mass  of  section  •  relative  to  the  vehicle  axes  -- 
arbitrary  values  given  prior  to  dynaaic  balancing. 

TAJTI  -  TAJIKS 
SEJKT  -  STRUCTURE 


products  of  inertia  of  section  i  referred  to  sectional  axes, 
as  in  ’/PCS. 

TAXR7S 

soda1,  inertia  terns. 

SFLSTS 


components  of  orthogonal  unlt_vector3  giving  directions  of 
acceleration  oriented  axes,  X  end. -Tit  being  parallel  to  the 
surface  of  the  fuel  in  tank  i,  and  1  being  perpendicular  to  the 
surface  of  the  fuel.  Same  os  in  VPCS. 

TALSRT 


,  _ 

L.  r  modal  inertic  terms  for  fuel. 
K  TALFTS 


modal  inertia  term3  for  structure. 


fM  p:  totd  masses  of  fuel  in  tanks.  Same  as  in  ’/PCS. 
TAMF7S 


*  * 


-c-  • 


•  v.'W- 


*  V  <  V 


t  .  -■/. 

-  •  •'*  ”’ 

*  i, 

<*  +■  *  «.  -  * V  •  ^  «» 


'3  o“  one  ir.e.*.js  tensor 


bending  moments  at.  -i  specified  1  c  o  a  t 
without  . i::d  (conpor.err.:  in  the  y  •: 
SPJrTTT 


;r.  flexible  ve.. 


bending  roncot  o.  flexible  ver. .  cle 
sF’i?Tr 


bending  r.onsr.t.  on  ri;;i 
force . . 

SPIRIT 


i  vehicle  ./itr.out  wind  and  without 


bending  moments  cn  rigid  vehicle  vitr.out  rind  but  with  thrust 
forces. 

SPMR2T 


bending  moments  cn  rigid  vehicle  with  .in.d  ano  with  thrus 
forces . 

SPHRir? 


aerodynamic  bending  moments  about  the  origin  Ju  to  elastic 
'w  deformation,  vithaut  *ind. 


SAMAET 


MA^  aerodynamic  bending  moments  about  the  origin  due  to  rigid 
w  motion,  -without  wind. 


SAMART 


M  A  m  aerodynamic  bending  moments  about  the  origin  due  to  elastic 
*•'  £ir“  deformation,  with  wind. 

SAKSBT 


A  aerodynamic  bending  moments  about  the  origin  due  to  rigid 

,vLn  RA-®  motion,  with  wind. 


SAMRBT 


\  J 


*  *  *  > 

-  . 


f\ 


.»  #■  *  *  ‘  *-  »  ‘  -  <  **  V 


*»  . 


■  ..  -  >  .) 

•  (  rf^f"  .  '  *'  x  •  •  /; 

♦  1  £‘  »  *  '  ‘  ■'  '  '  •*  >  •: 

'  ‘  r  -  •  -  ‘  •  •->  ■  •  ' 

V  ■  .  ♦  .  ;  \  V: 


-.1  ; 


M 


bonding  acnent;  about  the  or.  gin  due  to  gravity. 
ST-MCHT 


Wx  j.  inertial  bending  moments  about  the  origin  due  to  ela-tis  dor': 

"  nation . 

SIMI2T 


fvi  I  r  inert 

~  R®  SIM  IK? 


MTV  bending  moment:  about  the  origin  due  to  the  thrust  force:  of  the 
—  R«r  engine:. 


m  total  nass  of  vehicle  and  fuel  at  any  inotnnt.  Jane  as  in  VP  Cl. 

AMASS 


rn  i  nass  of  struct  iral  section  I-  Same  as  in  7FCS. 

TAMX7S 


m.\h  nass  of  the  H -til  particle  of  section  t. 
SEMIH 


effective  fuel  slosh  Basses  in  tank  !. 


wacs  of  part  of  the  vehicle  and  fuel. 
S1MP7S 


aass  of  port  of  section  i. 
SIMPIS 


**<*  »#*■»«-<  vo  s«v* 


phericai  and  lorizortil 


rr.  .  feel  -l.;'.  inertia  coupling  teras  for  5 
.*>*’. ir.dricnl  tank'. 

TAiffKS 

I 

t%!-  r.  "ter  of  “crodynaaic  part..  (cr  surfaces  }  in  section  i. 


-oral  inertia  tcr-r . 


N  sect! :nal  aerodyn'tuic  fcc-ading  •'.erse;  ternc  resulting  froa  elastic 

El  Set  motion,  •■•ithout  ;ind. 

5Al.r2XP 


(\j  .  3ectic-r.nl  aero iyr.a: .ic  bending  nonent  tei  •  resulting  froa  rigid 

fi  l  - .  .  ..  .S. - .  .J  .  1 

SANRIP 


sectional  acrcdyr.aalc  banding  coaent 
deformation,  'itr.  wine. 

SANHEP 


reiuitir  '  from  elastic 


N 


r 

RJrl 


.■ectionoi  aeroayra..ic  bending  :cuer.t  ttrr. 
notion,  ;ith  vind. 

SkFP.BT 


.-c-ilti.v 


ririi 


N  F  K  co-e_  ir.artio  torr.j  for  fuel. 

TAHFK5 

N  5  ncdal  inertia  terns  for  structure. 
*  SEEKS 


n 


nurfoer  of  elastic  degrees  of  free  ion. 
K02T. 
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/  f  ' 

cc’.pcr.^r.ts  in  the"V  cocidin.ite  system  cf  a  uni*  vector  at  osi.  t 
•A  Jhon  the  surface  of  sccticr.  .  ,  perpendicular  to  i.*.a  surface  end 
pointing  cutuard . 

AEr.HR 

P;  -.umber  of  r.s>*idcs  (or  ‘rs'e-eo)  ir.  acctitr  1 

NOP  I 


Ppr-t  products  jf  inertia  or  fuel  in  tar.i:  1  rc-ferv-fd  tc  axes  parallel 
to  vehicle  axes. 


Prij  p.odal  moments  and  negatives  of  products  of  lnertir  of  vehicle 
end  fuel . 

SEPJ7S  -  STRUCTURE 
TAPJ7S  -  TANKS 


r* 


sectional  coordinates  of  the  point  of  rotation  of  movable 
section.'.  Same  a..  in  7FC3. 

SEPPR 


R\ 


:::  liu.-  of  t?ir-;  I  •  n:s  it.  '/r»». 

TA^I'T 


radio.-  _f  equivalent  horizontal  cylindrical  tank  i  • 
Same  as  In  VPCS. 

TARirrr 


Rv;  radium  of  equivalent  vertical  cylindr< '—.1  t«oK  \  . 
Same  as  ^n  VPCS. 

TAKV7T 


\'i *■ 


nodal  p;  od'.’.ct'-  of  inertia  of  port  of  the  vehicle  referred  to 

/ehielo  axes. 

filRKRS 


.•aetisnnl  acvo-lyraut::  force  tor:..:  for  rigid  vehicle. 
AERTOT 


R'. 


r?.t> 

i 


r:  eat  tonal  nerc dynamic  -hear  force  ter:..;  for  rigid  vehicle. 
3 ARP IT 


% 

ir 

m 


& 


■1 


u  r*t>u- 


cectior.’.l  tare  Ayr  stile  handing  r.or.cr.i  tcrr.u;  icr  rigid  vehicle. 
CARPIT,  3ARP2T,  S  ARP  AT 


R 


r 

K 


SHIRKS 


i  c'-fe ;t  ratio-;  vf  "equivalent"  tanke.  Samo  as  in  VPCS. 

TARTL 

r a  quotient  cf  aspect  rutioa  of  rectangular  tank  i . 

TAOTL 
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!i t* wi  -fit  *r»*st 


s 


& 


’  rs 


v  •''.!. !o  ,.;.-u’  V.'-.i  •• 


i:.jrMr.  ..f 
i:  ,ri’c:;. 


tructure  of  -ai.fi -  j 


S»J 


* :  r  -  J  Cn 
AXSK7P 


i  -tn  -M. 


$1 

5  i 


c>r  chor.r  forcr;  .V  :i  j^ov or*  flexible  vehicle  uitheut 

,1M  (co.upor.unr  lr  ‘he  .•  r.cr:i*.;tc  c;v* -./.  '  & 

3PST7P  ’’ 


2  r  shear  forces 

*r  5P537P 


(•**..  i  «  a  -  ,  .Uf  .  (U.  .  »  „  ' 


li 


•hear  fovea:  due  ‘o  das4  i  :  i. 
SP527P 


•:*.>r;.«vior.  iltl.au*  find. 


I 

,ri 

:i 


—  i 

O  |  shear  force  3  due  tc 

e spr.SDP 


o'  .  tie  ddorr.'.atism  ith  ,/ind. 


shear  force:  or.  ri™12  vchi  .la  without  :in-l  and  ./ithout  thrust  force. 
R‘  SPSIllP 


R2 


shear  forces  on  rijitl  •.'chicle  .ithout  \r ind  but  oith  thru3t  forces. 
SPSR2P 


S  „  .  Jhaur  force::  or.  riciJ  vehicle  fith  ••■ind. 

R "  SPSRBP 
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S  3  'cctionsl  aerc Jyiiamlc  st-'ffness  terns. 

J*-v  A'-JSTuT 


aerodynamic  shear  forces  due  to  elastic-  jeforaation,  without  wind. 
~  StSRAP 


aerodynariic  sr.oor  forces  iue  to  rigid  notion,  without 


S^A£jr  aerodynamic  shear  forces  due  to  elastic  deformation,  with  wind. 


r- 

S  A  .  eerodynanlc  3hear  forces  due  to  rigid  notion,  ulth  wind. 
w  * v  3ASRBP 


S  A  El  elastic  contribution  of  section  i  to  the  ►erolynmoic  shorr 
^  forces,  without  wind. 

SASAHP 


S?A  .  rigid  contribution  of  section  i  to  the  aerodynamic  s,hear  forces, 
without  wind. 

SASAR? 


S  A,-*  cl  a- tic  contribution  of  section  i  tc  the  aerodynamic  snear  forces, 
f -*r,with  wind. 

SASEHf 


&  A  .  .  "Igld  contribution  of  section  i  to  the  aerodynamic  shear  forces 
R*r,‘vith  wind. 

3ASRE? 


5/*  1 

vi  shear  forces  due  to  gravity. 
SGSGRP 


■  *-<*♦ 


51  .  inertial  shear  forces  due  to  elastic  deforaatli.-- . 
~  E  SISIEP 


n.  r 


inertial  shear  fujees  due  to  rigid  notion. 
3ISIRP 


SJ  shear  forces  due  to  the  thrust  forces  oi  the  eng'res. 
ST3TRP 


T  nuaber  of  tanks.  Sane  as  in  ’/PCS. 

N&’AH 


~p  conponents  in  the  y  coordinate  systea  of  the  thrv.st  force  at  the 

'  x.!.  i-tit  nozzle. 

2fTXZP(l) 


1“  .  components  in  the  y  coordinate  systen  of  the  thrust  force  at  the 
'i'  i-th  nozzle. 

2-fTXZP(2) 


T,  •,  conponents  la  the  y  coordiaat-  system  of  the  thrust.  force  at  the 
*jfy  i  -th  nozzle. 

2fTX2P(3) 


— _  T 

|  .  sectional  aerodynonic  daaplng  terns. 

J*1  ASTIJT 


T .  .  sectional  aerodynaalc  shear  force  terns. 
.SAT? IT 


T"  .  sectional  aerodynaalc  bending  noaeot  terns. 
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.  ♦  ‘  ‘  ‘out 


ire. 


'‘ercxiyr.ani-:  terns . 
**  ‘  -i  V?3T 


U  sectional  eerodynanic'  sheer  force  terao. 

«■'•  34UPIT 


"*-  s.t.u 

Vj  sectional  aerodynamic  bending  moment  terms. 

*•  SAUP1T,  SAUP2T,  SAUPB? 


Vi.  -c  cosine  of  angle  between  resultant  acceleration  and  axis  of 
cylindrical  tank  i.  Sane  os  ii,  VPCS. 

TAUS 


\/  r  components  in  the  y  coordinate  system  of  the  linear  velocity  of 
the  vehicle  at  the  origin  of  the  vehicle  axes. 

A3VR7T 


*  1 

^ 


c opponents  of  the  velocity  of  the  wind. 
A3VRAT 


r  r  .  r 

v;  v-v‘- 

asvrst 


V  components  of  the  vehicle  velocity  at  the  center  of  mass. 

C  ASVBCf 

v;  dv‘A 

srvDcr 

i 

j/\j-  *  components  of  the  velocity  of  particle  A  of  section  \  relative  to 
the  vehicle  axes. 


ADY5KT 


VV  inertia  coupling  terns  for  part  of  the  veh4.-  le 

*  -  .,xes. 


:  j  „  ».-.n  tern  of  W 

<N  V. 


W«.  . 


UU'-(  number  of  fuel  slosh  modfs  in  each  direction  for  tenk  i  .  (.  —  ^  ) 
Kjfal 


AAJ;*.  the  "piston  speed"  (or  aovcwash)  at  the  J»-th  surface  of  the  i  -th 
section . 

AEV-THT 


ULT_r  o  oris  tern  cf  A  i;j, 
ABVSltT 


distance  along  the  axis  censured  fron  the  center  of  cylindrical 
tank  i  ,  nondlaenalonalized^  vlth  respect  to  the  tank  lee^-h.  T" 

U3e  vlth  the  real  tank,  'X.I  is  positive  la  the  direction  so  chosen 
In  the  WCS  data  to,  be  subaitted,  cunber  5*  In  use  vlth  a  vertical 
equivalent  tank,  Xtis  positive  in  the  direction  of  the  "resultant 
acceleration*  if  ui  la  positive,  end  positive  in  the  direction 
opposite  to  the  "reaultaol  acceleration"  ir  u;  li  negative..  (  ut « 
co*  Os  and  in  the  case  of  a  vertical  "equivalent”  tank  |u.l  >  I  /*%'.) 


X.  .  coordinates  of  geometric  center  of  tsnk  i  or  of "the  point  of 
'  >o Lotion  of  coracle  section  i  .  naae  as  in  VFCS. 


V<. 


Xr 

K  K\ 


*  * 

£  *  Hi  . 

*•1 

7JUSX  -  TAMS 
32XRK  -  STRUCTURE 


A  _  . 


Vi* : 


_»  r. 


-  p/r  =  VI  *i.  z  c 

r  it*  i  *  * 


one  tens  of  A  tii»  • 
ADXSH? 


static  unbalances  of  part  o'  the  vehicle  ref'’  -rod  to  vehicle 

axes. 

sr/RTs 


static  unbalances  of  part  of  section  t  ,  referred  to  vehicle 
axes. 

SIYRIS 


dynarslc  unbalances  of  part  of  the  vehicle  referred  to  vehicle  axes. 
SIYRK3 


inertia  coupling  terns  for  part  of  the  vehicle  referred  to 
vehicle  axes. 

SXYXXS 


certain  suaeatlons  connected  with  fuel  slosh  in  vertical 
cylindrical  tanka. 

TATPI-  TA^PP 


coordinates  of  the  h  -tb  particle  of  the  »  -th  section  in 
the  i jr  coordinate  aystea. 

AOfTifr,  SATKET 


products  of  inertia  of  part  of  section  *  ,  referred  partly  to 
sectional  axes  and  partly  to  vehicle  axes. 

srms 


aodal  products  of  inertia  of  part  of  section  i  ,  referred  partly 
to  sectional  axes  and  partly  to  vehicle  exes. 


.  (  - 


>  *'* 


•  vv  - d 


WSM*- 


.... 


7 

v. »  4— ‘-At>‘ 


tv 

“V' 


&KI 


distance  in  tank  *  fron  fuel  center  '■»»  buss  to  spring 
fflas  >  s ,  positive  up. 

TAZX7T 


distance  fron  gectnetric  center  to  center  of  fuel  mass, 
positive  upward,  for  equivalent  .ank  5.  Sane  at  in  VPC 
TAiBAP. 


coordinates  of  the  center  of  ss 
^CS. 

SB2SI 


:-p  of  settlor  i  .  Sane  as  in 


-5,  coordinates  of  tenter  of  mass  of  fuel  in  tans.  ».  3a*e  30  in 

A  Fi  '/PCS  - 

V  TA2FRT 


coordinates  of  the  centerof  mass  of  the  vehicle.  Same  as  in 
YPCS. 

ZCFHT 


cofsponents  of  partial  angulur  velocity  with  respect  to 

of  the  J  ft  coordinates  relative  to  the  J system  -  value* 

obtained  after  dynamic  balancing- 

SFAFR 

TAAPIl 


fuel  height  angle  for  spherical  and  horizontal  cylindrical  tanks; 
that  is,  the  angle  between  the  free  surface  and  a  lino  fron  the 
center  of  the  t^nk  to  the  intersection  of  the  free  surface  with 
the  wall  of  the  tank.  Same  as  in  VPCS. 

TA3KT3 


components  in  the  J *  syrten_^f  the  partial  angular  velocity _ 
with  respect  to  q*  of  the  «)  ri  coordinates;  relative  to  the  Jr 
system  —  arbitrary  values  given  prior  to  dynamic  balancing. 
SEBKI  -  STRUCTURE 
TA3KI  -  TANKS 


i* *.?  -  —  -  ♦v.  v *  £  -V" • 


.^‘u. **  ‘  Iff- 


«»  » 


\  •  ;  • 


uO'.ponents  in  the  J  r  system  of  the  partial  euag’.J.ar  velocity 
with  respect  to  qJ  of  the  coordinates  relative  to  roe 

Jr  aysteo  arbitrary  values  given  prior  to  dynamic  balancing. 

IABRK  -  TAIiiCl 
3SSPK  -  STHVC?UHOS 

products  of  inertia  of  section  i  ,  ref»  ere*  ,o  !he  sectioral 
axes . 

taccpo  -  ta:;k.j 


products  of  inertia  of  part  of  socMt.  i  ,  t-.red  to  the 
sectional  axes. 

3 ICO IS 


Y  Y  y  constants  obtained  from  Bessel  functions  and  used  with  vertlsax 
"  i  ,  "a  .  '•3  cylindrical  tanks. 


A  F 

JK 


inertia  coupling  terms. 


inertia  cotmllna  tarmc  for  fhel . 
7ADFJS 


AS..  inertia  coupling  terns  for  structure. 

'  33D5J3 


tine  increment  used  in  the  numerical  integration. 


nodal  inertia  terms. 
SPSPSS 


inertia  counl  ing  terms. 
O'^TTS 


inertia  coupling  terns. 


,*t>‘  -  -•  *  ~  f-k.  ‘  .* 


1 

f  iv_  __ 


\ 


9~Fok 


inertia  coupling  t'-rras  Tor  fuel. 
TATMFS 


inertia  coupling  terms  for  struct 
SKTASS  '  " 


m/vio  1 

TATVJO 


e-5o 


modal  inertia  term:;  for  structure 
SET3JS 


modal  inertia  terns. 


./TS* 

A., 


modal  products  of  inertia  of  section  ». 
7ACLAM  -  TAJIKS 
S2CMM  -  siaucruR5; 


modal  products  of  inertia  of  pert  of  section  •• 
3ICLL5 


AF 

AS 


modal  inertia  terms  for  fuel. 
TACtfS 


modal  inertia  terms  for  structure. 


^  frequency  parameters  for  lateral  motion  of  horizontal 
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